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ABSTRACT 


This  study  was  undertaken  to  determine  if  the  US  Army's  Ml  19 
Howitzer  could  be  safely  towed  in  the  “firing”  position  using  the 
existing  A-Frame  brace  to  support  the  gun  tube  during  movement.  A 
computer  model  of  the  howitzer  was  developed  that  would  predict 
motions  of  the  howitzer  and  consequently,  loads  on  critical 
components.  The  results  obtained  from  this  computer  model  were 
validated  by  comparing  them  to  actual  data  obtained  from  an 
instrumented  test  done  on  the  howitzer  at  Aberdeen  Proving 
Grounds,  Md.  The  validated  model  was  used  to  predict  the  "worst 
case"  loads  on  the  howitzer’s  A-Frame  and  the  T-Bar  during 
movement  in  the  "firing"  position.  The  maximum  predicted  stresses 
in  the  A-Frame  and  the  T-Bar  from  the  "worst  case"  loads  were 
compared  to  the  yield  stress.  These  comparisons  showed  that  the  A- 
Frame  and  T-Bar  were  strong  enough  to  withstand  the  "worst  case" 
predicted  loads  that  would  be  caused  by  towing  the  Ml  19  Howitzer 
in  the  "firing"  position. 
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I.  INTRODUCTION 


A.  GENERAL 

The  purpose  of  this  study  was  to  determine  the  feasibility  of 
transporting  the  Ml  19  howitzer  in  the  firing  position  over  terrain 
and  at  speeds  used  during  normal  tactical  operations.  The  end 
objective  was  to  determine  if  costly  design  changes  would  be 
required  to  safely  permit  towing  the  Ml  19  Howitzer  in  the  firing 
position.  The  approach  used  was  to  model  the  vertical  plane 

dynamics  of  the  towing  vehicle  (HUMMV)  and  the  howitzer,  model 
representative  terrain  inputs  to  the  wheels  in  the  form  of  potholes, 
and  to  perform  experiments  to  validate  the  model.  It  is  proven  that 
the  simulation  model  is  reasonable  and  it  follows  that  stresses  in  the 
critical  components  are  not  expected  to  warrant  redesign. 

B.  BACKGROUND 

The  Ml  19  Howitzer  is  a  British  designed  howitzer  that  is 
currently  being  built  at  Rock  Island  Arsenal  (RIA),  Illinois  and 
fielded  to  US  Army  Field  Artillery  units.  The  Ml  19  Howitzer 
operational  manuals  specify  towing  procedures  that  require  the  gun 
tube  to  be  placed  in  the  stowed  position  when  moving  the  howitzer 
during  tactical  operations.  These  towing  procedures  are  different 
and  more  difficult  to  implement  than  those  used  for  all  other  towed 
howitzers  in  the  US  Army.  All  other  US  Army  towed  howitzers  are 
towed  during  tactical  operations  with  the  gun  tube  in  the  firing 
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position,  facilitating  rapid  emplacements  and  displacements  that  are 
required  for  successful  tactical  operations. 

The  Ml  19  Howitzer  has  an  existing  A-Frame  brace  which  can 
support  the  gun  tube  when  the  howitzer  is  towed  in  the  firing 
position.  The  primary  concern  of  this  study  was  whether  or  not  the 
A-Frame  and  the  T-Bar  structure  on  the  howitzer’s  trail  section,  to 
which  the  bottom  of  the  A-Frame  brace  connects,  will  be  strong 
enough  to  support  the  loads  generated  when  towing  the  howitzer 
over  rough  terrain.  This  study  also  addresses  concerns  about  trail 
loads  on  the  howitzer.  This  study  does  not  evaluate  rollover  stability 
for  the  Ml  19  Howitzer  when  towing  in  the  firing  position. 

The  terrain  that  the  HUMMV/Howitzer  system  can  traverse  is 
clearly  limited  by  the  terrain  the  HUMMV  can  negotiate,  therefore, 
the  Ml  19  Howitzer  and  HUMMV  have  been  evaluated  as  a  system.  A 
dynamic  simulation  model  of  this  system  was  seen  to  be  a  necessary 
tool  to  aid  in  making  design  decisions  from  a  rational  viewpoint  and 
is  the  basis  of  this  study. 

C.  METHOD  OF  ANALYSIS 

To  analyze  the  Ml  19  Howitzer  a  computer  model  of  the 
HUMMV/Howitzer  system  was  developed  to  predict  the  response  of 
the  system  to  roadway  inputs.  This  computer  model  treats  the  major 
mass  elements  of  the  HUMMV/Howitzer  system  as  lumped  masses  to 
develop  vertical  plane  heave  and  pitch  equations  of  motion  for  the 
system.  The  development  of  this  computer  model  is  contained  in 
Chapter  II.  To  drive  the  computer  simulation  model,  terrain  models 


2 


were  developed  for  both  a  “worst  case"  trench  type  pothole  and  a 
sinusoidal  "washboard"  roadway.  Details  on  the  development  of 
these  roadway  inputs  are  contained  in  Chapter  III.  Chapter  IV 
contains  the  results  of  driving  the  computer  model  over  the  roadway 
inputs.  Chapter  V  explains  the  structural  analysis  that  was  used  to 
determine  the  maximum  stress  in  the  A-Frame  and  the  T-Bar. 
Chapter  V  also  contains  the  procedures  and  results  of  calibrating  an 
A-Frame  that  was  used  as  a  load  cell  for  an  instrumented  dynamic 
test.  Chapter  VI  explains  the  conduct  and  the  results  of  the 
instrumented  test  on  the  Ml  19  Howitzer/HUMMV  system  that  was 
conducted  at  Aberdeen  Proving  Grounds,  Maryland  in  April  1993. 
Also,  in  Chapter  VI  the  results  of  the  model  are  compared  to  the 
experimental  results  from  the  instrumented  test.  This  comparison 
showed  reasonable  correlation  between  the  model  and  experimental 
results.  Conclusions  and  recommendations  based  on  this  study  are 
contained  Chapter  VII.  The  final  recommendation  is  that  towing  in 
the  firing  position  should  be  allowed. 
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11.  nUMMV/HOWITZER  SYSTEM  COMPUTER  MODEL 


The  computer  model  developed  is  based  only  on  the  heave  and 
pitch  equations  of  motion  (EOM)  for  the  HUMMV/Howitzer  system. 
The  reason  for  this  simpliflcatioo  is  that  the  lateral,  yaw,  and  roll 

motions  are  much  less  significant  for  determining  the  forces  in  the 
A-Frame  for  the  type  terrain  inputs  considered  by  this  study.  While 
horizontal  plane  motions  are  key  elements  for  the  analysis  of  turns 

and  rolling  during  transit,  the  concerns  of  this  study  indicate  that  a 
vertical  plane  model  is  adequate.  The  major  mass  groups  of  the 
HUMMV/Howitzer  were  treated  as  lumped  masses  in  the 
development  of  the  equations  of  motion.  The  equations  of  motion 
were  developed  using  Lagrange  techniques. 

A.  HUMMV/HOWITZER  SYSTEM  PARAMETERS 

The  primary  sources  of  information  for  system  parameters  for 
the  HUMMV  and  howitzer  were  the  US  Army  Tank- Automotive 
Command  Research,  Development  &  Engineering  Center  Technical 
Report  No.  13337,  [Ref.  1],  for  HUMMV  data  and  US  Army  Combat 
Systems  Test  Activity,  APG,  Report  No.  P-83179,  [Ref.  2],  for  the 
Ml  19  Howitzer  data.  Additionally,  many  system  parameters  were 

determined  from  direct  measurement  on  the  HUMMV  and  howitzer 
or  from  design  diagrams,  or  obtained  from  the  study  sponsor  at  Rock 
Island  Arsenal  (RIA).  For  readers  of  this  study  that  are  unfamiliar 

with  the  components  of  the  Ml  19  Howitzer,  Appendix  G  contains 
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several  diagrams  that  depict  the  howitzer  components  that  are 
important  for  this  study. 

The  system  parameters  used  for  the  HUMMV  are  based  on  using 
the  Model  No.  M1069  HUMMV  as  the  prime  mover  for  the  Ml  19 
Howitzer.  The  cargo  loading  of  the  HUMMV  signiHcantly  effects  its 
center  of  gravity  and  consequently  effects  all  of  the  HUMMV's 
system  parameters  which  are  dependent  on  its  center  of  gravity. 
Modeling  was  therefore  limited  to  using  the  same  cargo  loading 
configuration  that  was  used  in  [Ref.  1]  and  other  loading  effects  were 
ignored.  These  HUMMV  simplifications  are  justified  because  for  this 
computer  model  the  HUMMV  serves  primarily  as  the  driving  force  at 
the  howitzer  hitch  causing  pitching  motions  of  the  entire  howitzer 
and  we  are  not  concerned  with  the  loading  of  HUMMV  components. 
The  major  components  of  the  HUMMV/Howitzer  system  computer 
model  are  identified  in  Figure  2-1  below. 

In  the  model,  the  tires  act  as  linear  springs  which  provide  force 
inputs  to  the  HUMMV  and  howitzer  wheels,  based  on  the  terrain  over 
which  the  model  was  "driven".  The  suspension  systems  for  both  the 
HUMMV  and  the  howitzer  are  modeled  as  springs  and  dampers 
connected  in  parallel.  The  A-Frame  and  T-Bar  are  modeled  as  a  set 
of  stiff  springs  connected  in  series.  The  following  paragraphs 
provide  more  detail  on  the  modeling  of  the  individual  components 
and  the  determination  of  appropriate  parameter  values.  Appendix  A 
provides  a  listing  of  the  numerical  values,  units,  and  the  symbols 
used  for  each  system  parameter  in  the  computer  model.  The 
parameter  values  for  those  components  with  identical  components  on 
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the  right  and  left  sides  of  the  system  (wheels  and  suspension 
components)  are  multiplied  by  two.  The  values  in  Appendix  A  for 
these  components  reflect  the  effects  of  both  sides. 

1.  HUMMV  and  Howitzer  Tires 

The  HUMMV  has  the  same  tires  on  the  front  and  rear, 
however,  the  front  and  rear  tire  pressure  is  different  (20  psi  vs  30 
psi)  resulting  in  the  rear  tire  spring  being  much  stiffer.  The 

stiffness  used  for  the  HUMMV  tires  was  obtained  from  normal  force 
vs  tire  deflection  data  [Ref.  1].  The  stiffness  of  the  howitzer  tire  used 
in  the  model  was  estimated  because  no  measured  data  on  the 
howitzer  tire  stiffness  was  available. 

The  tires  are  modeled  as  linear  springs  which  transmit  a 
vertical  force  to  the  respective  wheel  masses  proportional  to  the 
amount  that  the  springs  are  stretched  or  compressed.  The  amount  of 
force  is  determined  by  the  relative  dynamic  displacement  between 
the  terrain  input  for  the  wheel  and  the  center  of  mass  of  the  wheel. 
If  the  tire  leaves  contact  with  the  ground  there  will  be  no  input  from 
the  terrain  until  it  again  returns  to  contact.  To  account  for  this  loss 
of  contact  situation  in  the  computer  model,  the  tire  spring  force  is  set 
to  the  static  weight  that  the  tire  supports  when  the  tensile  (de¬ 
compressive)  spring  force  exceeds  the  static  weight  on  the  tire.  The 
tire  springs  are  turned  back  on  in  the  computer  model  when  the 
conditions  are  met  for  road  surface  contact  to  be  re-established. 

2.  HUMMV  Suspension  System 

The  HUMMV  suspension  system  consists  of  a  shock 
absorber  and  coil  spring  for  each  wheel.  The  shock  absorber  passes 
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through  the  middle  of  the  coil  spring.  The  force  in  the  spring  is 

linearly  proportional  to  the  relative  displacement  between  the  center 
of  mass  of  the  wheel  and  the  attachment  point  on  the  HUMMV  body. 
The  values  for  the  spring  constants  were  obtained  from  Table  5-13 
of  [Ref.  1].  Precise  modeling  of  the  shock  absorbers  is  difficult 

because  they  have  a  mid-stroke  region  and  a  hydraulic  bump  stop 

region.  In  the  mid-stroke  region  the  shock  absorber  force  is  linearly 
dependent  on  the  relative  velocity  of  the  center  of  mass  of  the  wheel 
and  the  attachment  point  on  the  HUMMV.  However,  the  linear 

damping  constant  changes  depending  on  whether  the  force  in  the 
shock  absorber  is  tensile  or  compressive.  When  the  shock  absorber 
is  compressed  or  stretched  into  the  hydraulic  bump  stop  region  the 

damping  of  the  shock  absorber  greatly  increases.  This  sudden 

increase  in  damping  causes  problems  in  the  integration  routine  that 
was  used  for  the  computer  model.  To  overcome  these  problems  and 

still  account  for  the  energy  loss  due  to  the  increased  damping  in  the 
hydraulic  bump  stop  region,  the  mid-stroke  linear  damping  constant 
obtained  from  [Ref.  1]  was  increased  by  20%  for  both  the  front  and 

rear  shock  HUMMV  absorbers.  Figures  2-2  and  2-3  shows  the  force 
versus  velocity  plots  that  were  obtained  from  [Ref.  1]  for  the  front 
and  rear  shock  absorbers,  respectively. 

3.  Howitzer  Suspension  System 

Modeling  the  howitzer  suspension  system  is  more 

complicated  than  the  HUMMV  because  it  consists  of  a  torsional  spring 
and  a  shock  absorber  that  is  not  vertically  oriented.  The  torsional 
spring  is  modeled  as  a  linear  vertical  spring  that  is  connected  to  the 
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center  of  mass  of  the  howitzer  wheel  and  attached  to  the  howitzer 

trail  directly  above  the  wheel  center  of  mass.  The  effective  spring 

stiffness  constant  was  calculated  based  on  the  measured  data  of  the 
torsional  spring  stiffness.  Figures  2-4  through  2-7  show  the  angular 
deflection  versus  vertical  force  plots  that  were  contained  in  [Ref.  2]. 
Using  the  average  torsional  spring  stiffness  from  this  data  and  the 
length  of  the  road  arm,  an  angular  rotation  was  converted  to  an 
equivalent  vertical  displacement  at  the  end  of  the  road  arm.  The 

spring  constant  for  the  model  was  then  calculated  by  determining 

the  constant  value  that  must  be  multiplied  times  this  vertical 

displacement  at  the  end  of  the  road  arm,  to  produce  a  moment  that  is 

equivalent  to  that  caused  by  the  torsional  spring  for  a  specified 

angular  rotation. 

A  similar  procedure  was  used  to  model  the  howitzer's 

shock  absorber  as  a  vertical  shock  absorber  with  a  linear  damping 
constant  that  is  velocity  dependent.  The  attachment  points  for  the 

ends  of  the  shock  absorbers  in  the  model  are  the  same  as  the 

modeled  end  points  of  the  howitzer  suspension  spring.  The  actual 
non-vertical  shock  absorber  force  was  transformed  to  a  vertical  force 
in  the  model  by  determining  the  amount  of  vertical  force  applied  a  t 
the  end  of  the  road  arm  needed  to  produce  an  equivalent  moment 

about  the  road  arm  pivot  point.  Like  the  HUMMV  shock  absorbers, 
the  howitzer  shock  absorbers  are  non-linear  and  have  different 

damping  values  for  tension  and  compression  that  are  accounted  for 
in  the  computer  model.  This  non-linear  damping  in  the  shock 
absorbers  is  a  significant  aspect  of  the  computer  model.  The  values 
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used  for  the  shock  absorber  damping  constants  were  determined 
from  force  versus  velocity  test  data  provided  by  RIA.  Figure  2-8 
shows  the  force  versus  velocity  plot  for  the  howitzer  shock  absorber. 


Ml  19  HOWITZER,  105mm  GUN,  Vehicle  No.  27 


Figure  2-6  Torsional  Spring  Stiffness  for  Serial  No.  27  (Upward 


Motion) 


efiection  -  degrees 
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Figure  2-7  Torsional  Spring  Stiffness  for  Serial  No.  27  (Downward 
Motion) 
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Figure  2-8  Force  versus  Velocity  Plot  for  Howitzer  Shock 
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4.  A-Frame/T-Bar  Stiffness 

The  A-Frame  and  T-Bar  are  modeled  as  stiff  springs 
connected  in  series  to  determine  a  combined  linear  spring  stiffness. 
The  stiffness  of  the  A-Frame  was  determined  during  the  A-Frame 
calibration  test  (explained  in  Chapter  V).  The  stiffness  of  the  T-Bar 

was  determined  by  using  the  principle  of  superposition  of  forces  to 
determine  the  effective  bending  stiffness  for  the  non-symmetric 
bending  of  the  T-Bar  that  is  caused  by  the  A-Frame.  This  non- 

symmetric  bending  of  the  T-Bar  results  because  the  angle  at  which 
the  A-Frame  transmits  forces  and  moments  to  the  T-Bar  is  not 

perpendicular  to  an  axis  of  symmetry  of  the  T-Bar. 

5.  Equilibrator  Springs 

The  equilibrator  springs  provide  only  small  dynamic  forces 
because  their  combined  stiffness  is  nearly  1000  times  less  than  the 
effective  stiffness  of  the  A-Frame/T-Bar  assembly.  Therefore, 
equilibrator  springs  are  neglected  in  the  computer  model. 

6.  Elevating  Clutch 

The  effects  of  the  elevating  clutch  of  the  howitzer  were 

neglected  in  the  flnal  computer  model.  It  is  believed  backlash  in  the 
elevating  mechanism  gears  will  prevent  the  clutch  from  being 
activated  while  towing  with  the  A-Frame  attached  because  the 
relative  rotation  of  the  gun  tube  is  very  small.  However,  because  of 

some  uncertainty  about  the  clutch  several  computer  runs  were  made 
with  clutch  modeled  as  a  linear  torsional  damper.  The  value  of  the 
torsional  damper  used  in  the  model  was  estimated  because  no  test 
data  on  the  clutch  was  available.  These  model  runs  with  the  clutch 


included,  produced  only  slightly  different  results  from  the  runs 

without  the  clutch  included.  Therefore,  because  of  the  previously 

mentioned  belief  that  the  clutch  would  not  be  activated,  it  was 
neglected  in  the  final  model. 

7.  HUMMV  Component  Mass  Properties 

The  HUMMV  is  divided  into  three  lumped  masses,  the 
HUMMV  body  mass  and  two  wheel  masses.  The  front  and  back 
HUMMV  wheels  are  treated  as  having  identical  mass  values  in  the 

model.  Each  wheel  weight/mass  includes  the  weight  of  the  tire, 
wheel,  shocks,  suspension  springs,  and  suspension  support  arms.  The 

HUMMV  body  weight/mass  varies  dependent  upon  it's  cargo  load. 
The  location  of  the  center  of  gravity  also  can  vary  and  is  determined 

by  the  positioning  of  the  cargo.  Therefore,  the  values  used  in  the 

computer  model  for  the  HUMMV  mass,  the  mass  moment  of  inertia 

for  pitch,  and  lengths  that  are  based  dependent  on  the  center  of 

gravity,  are  estimated  values.  These  estimated  values  are  based  on 

the  same  loading  situation  as  was  used  in  [Ref.  1].  The  pitch  mass 
moment  of  inertia  for  the  HUMMV  is  about  the  center  of  gravity  of 
the  sprung  HUMMV  mass. 

8.  Howitzer  Component  Mass  Properties 

The  howitzer  is  divided  into  three  separate  lumped  masses 
in  the  model;  the  howitzer  wheel  masses,  the  elevating  mass,  and 
the  combined  trail/saddle  mass.  The  wheel  weight/mass  includes 
the  tire,  rim,  and  approximately  30%  of  the  road  arm  that  rests  on 

the  wheel  assembly.  The  elevating  mass  includes  the  gun  tube,  the 

cradle,  the  recoil  system  and  all  other  parts  that  move  the  when  the 


gun  is  elevated  about  the  trunion.  The  mass  moment  of  inertia  of  the 
gun  tube  for  rotation  about  the  trunion  was  estimated  using  the 
approximate  masses  of  all  the  elevating  components.  The  center  of 
gravity  location  for  the  elevating  mass  was  estimated  by  RIA.  The 
combined  trail/saddle  mass  includes  all  of  the  gun  mass  except  the 
elevating  mass  and  the  wheel  mass.  The  center  of  gravity  location 
for  the  trail/saddle  mass  and  the  pitching  mass  moment  of  inertia  for 
the  trail/saddle  mass  was  estimated  from  design  diagrams. 

9.  Length  Properties 

All  lengths  used  in  the  computer  model  are  identified  in 
Figure  2-9  below.  These  length  properties  were  determined  either 
by  direct  measurements  on  the  howitzer  or  by  scaling  from  design 
diagrams. 

B.  EQUATIONS  OF  MOTION  (EOM) 

In  this  section  the  heave  and  pitch  equations  of  motion  (EOM) 
for  the  model  are  developed  from  Lagrange's  equations  using  the 
HUMMV/Howitzer  system  parameters  that  were  previously 
explained.  Lagrange's  equations  require  the  determination  of  the 
system's  total  kinetic  and  potential  energy  and  the  power  dissipated 
by  the  dampers,  in  terms  of  a  set  of  generalized  coordinates. 

Throughout  this  development  of  the  EOM,  expressions  are 
developed  that  are  functions  of  generalized  positions  and  velocities, 
angular  displacements,  angular  velocities,  and  angular  accelerations. 
To  simplify  the  EOM  with  minimal  loss  of  representational  accuracy, 
these  expressions  are  linearized  using  small  angle  approximations. 
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Figure  2-9  Length  Properties  of  HUMMV/Howitzer  System 


The  distance  from  the  howitzer  trunion  to  the  CG  of  the  Gun  Tube 
The  distance  from  the  howitzer  tmnion  to  the  A-Frame  connection 


The  sine  of  a  small  angle  is  approximated  as  just  the  angle  (sin6=6) 
and  the  cosine  of  a  small  angle  is  approximated  as  equal  to  one 
(cos6*I).  These  approximations  are  valid  because  the  pitch  or 
angular  rotations  are  predicted  to  be  less  than  five  degrees.  Order  of 
magnitude  considerations  are  used  to  further  simplify  the  EOM. 
Because  the  angular  displacements  are  small  their  angular  velocities 

and  angular  accelerations  are  also  small.  Therefore,  any  terms  that 

contain  an  angular  displacement,  angular  velocity,  or  angular 

acceleration  multiplied  times  another  angular  displacement,  angular 
velocity,  or  angular  acceleration  are  neglected  because  these  terms 

will  be  very  small. 

The  general  form  of  Lagrange's  equations  used  is: 

£L1  aT  ^  au  ^  ap 

dtU^J  Sqi  ^  5qj  aq.  ' 

where: 

qi  =  generalized  coordinates 
T  =  total  kinetic  energy  of  the  system 

U  =  change  in  the  potential  energy  of  the  system  with  respect  to 
its  potential  energy  in  the  static-equilibrium  position 
F  =  the  power  dissipated  by  forces  in  dampers 
Q  =  the  input  force  acting  on  the  system. 

The  development  of  the  above  components  of  Lagrange's  equation 
and  the  assembly  of  these  components  into  the  system's  heave  and 
pitch  EOM  follows.  Due  to  the  length  of  this  development  the  reader 
is  provided  only  the  final  results  from  each  step.  The  intermediate 
multiplication's,  differentiation's,  and  linearization's  are  omitted. 


1.  Generalized  Coordinates  (qj) 

The  set  of  generalized  coordinates,  qj,  that  was  used  to 
uniquely  specify  the  motion  of  the  HUMMV/Howitzer  system  are 
depicted  on  Figure  2-10  and  are  defined  follows: 

*  Z1&Z2:  The  vertical  displacements  for  the  front  and  rear 

HUMMV  wheel  masses,  respectively. 

*  Z3:  The  vertical  displacement  of  the  center  of  gravity  of  the 

HUMMV  body. 

*  Z4:  The  vertical  displacements  of  the  center  of  gravity  of  the 

howitzer  wheel  mass. 

*  01 :  The  pitch  angle  of  the  HUMMV  body  about  its  center  of 

gravity. 

*  62:  The  pitch  angle  of  the  howitzer  trail/saddle  mass  .about  its 

center  of  gravity. 

*  63:  The  absolute  rotation  of  the  gun  tube. 

As  depicted  on  Figure  2-10,  upward  vertical  motions  and  clockwise 
rotational  motions  are  defined  as  positive.  The  equilibrium  position 
for  each  of  the  generalized  coordinates  is  the  static  equilibrium 
position  for  the  respective  components  that  are  correlated  to  each  of 
the  generalized  coordinates.  In  defining  positions  to  key  locations  in 
the  system  unit  vector  notation  of  iy  j,  and  k  will  be  used 
corresponding  to  the  x,  y,  and  z  axis  directions  of  the  Cartesian 
coordinate  reference  system  shown  on  Figure  2-10. 

The  kinetic  and  potential  energy  of  the  system  and  the 
damping  and  input  forces  are  all  computed  in  terms  of  these 
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Figure  2-10  Generalized  Coordinates  for  the  HUMMV/Howitzer 
System 
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U3 :  The  road  Input  for  the  howitzer  wheel 


generalized  coordinates  and  the  road  inputs  U],  U2,  and  U3.  To  carry 
out  the  required  computations,  position  vectors  to  key  locations  on 

the  HUMMV  and  howitzer  are  defined  in  terms  of  the  generalized 
coordinates.  The  position  vectors  to  these  key  locations  are  depicted 
on  Figure  2-11  and  defined  as  follows; 

*  ri  &  ri  =  the  front  and  rear  attachment  points,  respectively,  for 

the  suspension  springs  and  shock  absorbers  on  the  HUMMV. 

*  r3=  the  center  of  mass  of  the  trail/saddle  mass 

*  r4=  the  modeled  attachment  location  of  the  suspension  spring 
and  shock  absorber  of  the  howitzer 

*  r5=  the  location  where  the  A-Frame  attaches  to  the  T-Bar 

*  r6=  the  location  of  the  center  of  mass  of  the  gun  tube/elevating 

mass 

*  r7=  the  location  where  the  A-Frame  attaches  to  the  gun 
tube/elevating  mass 

These  key  locations  that  are  defined  above  can  be  expressed  as 

position  vectors  in  terms  of  the  above  defined  lengths  (Figure  2-9), 
generalized  coordinates  (Figure  2-10),  and  the  unit  vectors  i,  j,  and  k 
of  the  Cartesian  coordinate  reference  system  (Figure  2-10)  as  follows; 

*  r,  =(Z,  -  L,sin(0,))k  +  (L,cos(0i))i 

*  Tj  =(^  +  L2sin(0i))k  -  (L2Cos(0())i 

*  “(Z,  +  LjSin(0j)  +  LnSin(02))k  -  (LjCos(0i)  +  L,iCOs(02))i 

*  -(2^  +  L3sin(0,)  +  L,2sin(02))k  -  (L,cos(0,)  +  L,2cos(02))i 

*  r,  “(Zj  +  L,sin(0,)  +  (L,3+ L,*)sin(02))k  -  (L3COS(0,)  +  (L,3 +  L,j)cos(02))i 


Figure  2-11  Key  Locations  for  the  HUMMV/Howitzer  System 
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s  the  location  where  the  A-Frame  attaches  to  the  gun  tube/eievating  mass 


♦  =  (Zj  +  L3  sin(e,)  +  LisSinCe^)  +  LcGsin(0o  +  0j))k 

*  r,  =  (Z,  +  Lj sin(0,)  +  L,sSin(02)  +  L^sin(0o  +  0j))k 


Lj  cos(0,)  +L,sCOs(02)Y 

Lj  cos(0,)  +  Li6Cos(02)Y 
+  LaCOs(0o +0j)  / 


The  velocities  of  these  key  locations  can  be  determined  by 
taking  the  time  derivatives  of  the  position  vectors.  However,  only 
the  velocities  of  locations  where  there  are  velocity  dependent  actions 
occurring  in  Lagrange's  equations  need  to  be  determined.  The 
required  velocities  of  the  key  locations  on  the  HUMMV  and  howitzer 
are  as  follows: 

*  r,  “(Zj  -  L,cos(0i)0,)k  -  (LiSin(0,)0,)i 

*  1*3  =(Zj  +  L2cos(0i)0,)k  +  (L2sin(0,)0,)i 

r3=(2^  +  L3Cos(0,)0i  +  L,icos(02)02)k  +  (L3sin(0,)0,  +  Li,sin(02)02)i 

*T^=[z^  +  L3Cos(0,)0,  +  Li2cos(02)02)k  +  (L3Sin(0j)0,  +  L,2Sin(02)02)i 

'Z3  +  L3COs(0i)0i  +  L,4Cos(02)02Y  rL3sin(0,)0,  +  L,jSin(02)02\ 

*  r*  =  k  +  I 

\^  + l^r^G  C®s(0o  03)03  J  L(.oSin(0o  +  03)03  j 


In  addition  to  the  above  velocities,  the  vertical  velocity  of  the  CG  of 
the  HUMMV,  Z3,  and  each  of  the  wheel  masses,  Z,,  Z2,  &  Z^  are 

needed  in  Lagrange's  equation. 

2.  Kinetic  Energy  (T) 


The  velocities  defined  above  are  used  to  calculate  the 
kinetic  energy  of  the  system.  The  kinetic  energy  is  calculated  by 
summing  the  linear  and  rotational  kinetic  energy  for  each  of  the 


system's  masses.  The  total  kinetic  energy  of  the  system  is  given  by: 
£  T,  -  £  i  -  1,6;  j  -  1,3 

where: 

Ti  =  the  kinetic  energy  of  each  mass  in  the  system 
m  i  =  the  mass  of  each  lumped  mass 
Vj  =  the  linear  velocity  of  the  each  mass 
Ij=  the  heave  mass  moment  of  inertia  of  each  rotating  mass 
about  its  center  of  gravity 

0j=  the  angular  velocity  of  each  rotating  mass 

Therefore,  the  total  kinetic  energy  of  the  HUMMV/howitzer  system 

was  determined  by  summing  the  kinetic  energy  components  as 

follows: 

+  ^1,(6,)"  +  iij(e,)=  ♦  ii,(e,)’ 

In  the  above  expression  the  square  of  the  velocities  is  determined  by 

substituting  in  the  appropriate  velocity  expression  as  defined  above, 

and  taking  the  dot  product  of  the  velocity  with  itself.  This 

substitution  and  dot  product  operation  is  left  for  the  reader.  When 

taking  this  dot  product  all  of  the  k  dot  i  and  i  dot  k  cross  product 

terms  drop  out  because  the  i  and  k  unit  vectors  are  orthogonal. 

Having  determined  the  total  kinetic  energy,  the  time  derivative  of 

the  partial  derivative  of  the  total  kinetic  energy  with  respect  to  the 

Q  T 

velocity  of  each  of  the  generalized  coordinates,  — ■—  can  be 

dtdq, 

determined.  Carrying  out  this  differentiation  and  linearizing  yields 
the  following: 


^  d  dT 
dtaZi 

*1£L 

dtazj 

£  ar 
*  dt  az. 


M,Z. 

MjZj 

(M,+M5  +  M,)Z3  +  (M5L3  +  M,L3)e;  +  (MjL,,  + 


+  MjL(.gCOs(6Q)63 


*liL 

dtaz^ 
d  ar 

*  dt  a0j 


M,Z, 

(MsL3  +  MjL3)Z3  +  (I,  + 


M3L3^+  M,L3^)e. 


+  (MjL3Lm  +  MjL3L|^02  +  M^L3LpQ  cos(00)03 

J  arp 

*  (Jt  00  *  ^5^11  ^6^16)^  (^5^3^11  ■*' 

+  (I2  +  +  M0L,0LegCOS(0o)03 

d  ax 

“  M6LgQCOs(0Q)Zj  +  M0L3LgQ  cos(0o)6]  +  M0Lj0LgQCOS(0o)02 

*  dt  UO3 


+  (I3  +  M^Lco")©; 


The  total  kinetic  energy,  T,  must  also  be  differentiated 
with  respect  to  each  of  the  generalized  positional  coordinates, 
however,  for  this  system  this  differentiation  is  equal  to  zero  for  all 
coordinates. 

3.  Potential  Energy  (U) 

The  potential  energy  for  this  system  results  from  the 
strain  energy  that  is  developed  in  the  suspension  springs  and  the  A- 
Frame  spring.  The  gravitational  potential  energy  that  arises  from 
changes  in  vertical  positions  is  ignored  because  the  resulting  terms 
drop  out  of  the  dynamic  equations  of  motion.  The  potential  energy 
in  any  spring  is  determined  as  follows: 
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U,-iK,(Ax,)^ 


where: 

Ui  =  the  potential  strain  energy  in  each  spring 

Kj  =  the  spnng  stiffness  for  spring 

Ax.  =  the  change  in  length  of  the  spring. 


Since  locations  to  the  attachment  points  of  the  ends  of  all  the  springs 
have  been  defined  above  in  a  common  inertial  coordinate  reference 
system,  taking  the  difference  of  the  position  vectors  of  the  end  points 
of  each  spring  will  yield  the  change  in  length  of  the  spring.  Thus  for 
this  system,  the  total  potential  energy  is  found  by  summing  the 
potential  energy  in  each  of  the  springs  as  follows; 


iK3(Z2-Ujf  » 


Therefore,  the  partial  derivative  of  the  total  potential  energy  with 

respect  to  each  of  the  generalized  coordinates,  — ,  after  linearizing 

aq; 

results  in  the  following: 


*  —  -  K.Z.  +  K,Z.  -  K3Z3  +  K,L,8,  -  K,U. 

az, 

*  It  *  '‘>^2  *  ■  •^'^2 

aZj 

-K,Z,-  K,Zj  +  (K2  +  K,+K,)Z3  -  K,Z, 

*  aZj 

+  (KjL,  +  K4L2  +  K  74)01  +  K  74202 
-K,z,  +  (K.tK,)Z.  -  K,L,8,  -  K,L,JB, 

OZ^ 


K*U3 
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au 

*  aOi 


au 
*  ae^ 


-  KjLiZ,  -  K4L2Z2  +  (-K2L,  +  K,Lj  +  K7L3)Z3  -  K7L, 

(KjL.^  +  +  K,L,^)e,  +  K^LjL,^, 

"  KyLijZj  ~  K7L,2^4  ■*■  ■•■  KyLjj  ®2 

+  KALAL,3COs(eo)02  -  K^L^L, 3008(60)63 
*  3  005(60)62  +  K^Ly^Lj300S(6o)63 


4.  Generalized  Power  Lost  Due  to  Damping  Forces  (F) 

The  amount  of  generalized  power  that  is  dissipated  by 
eaoh  of  the  dampers  in  the  system  given  by  the  equation: 


F,  -  iBj(Ax,)" 

where: 

Fi  =  the  generalized  power  dissipated  in  eaoh  damper 

Bi  =  the  damping  constant  for  each  damper 

AXj  s  the  difference  in  velocity  of  the  attaching  points  of  each 

damper. 


The  total  amount  of  generalized  power  dissipated  in  the  system  is 
found  by  summing  the  energy  dissipated  in  each  dampers.  Summing 
the  energy  dissipated  by  the  damping  forces  yields  the  following: 

F-iB,(r,  -  Z,f  *  iB,(r.  -  Z.)’  t  -  Z,f  *  iB,(0,  -  e,f 


In  Lagrange's  equation  this  total  dissipated  energy,  F,  must  be 

differentiated  with  respect  to  the  velocities  of  each  of  the  generalized 
dF 

coordinates,  — r-. 

aqj 
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Carrying  out  this  differentiation  and  again  linearizing  results  in  the 
following: 


-  ^  -  B,Z3  +  B,L.e. 

*  B,Z,  -  -  B,L,e, 

*  -B,Z,  -  B,Z,  ♦  (B,  +  B, ♦  B,)Z,  -  B,Z.  t  (-B,L,  +  8,4)6, 

'  W' 

(®l^  *  64^)^  “  ®,1-S 

*  (/tJj 

+  (8,4=  t  B,L,“  ♦  B,L,')e,  »  B,L,4J9, 

*  M:  *  ^  B.)4  -  BA 

*  g-  -BA  .  BA 


5.  External  Forces  (Qi) 

The  external  forces  acting  on  the  system  are  the  roadway 
inputs  that  are  transmitted  through  the  tire  springs.  For  the 
computer  model  these  input  forces  are  treated  as  purely  vertical 
forces.  The  forces  result  from  the  road  inputs  causing  a  change  in 
length  of  the  tire  springs.  These  road  inputs  are  depicted  on  Figure 
2-10  and  defined  as  follows  : 

*  Ui  =  the  road  input  for  the  front  HUMMV  wheel 

*  U2  =  the  road  input  for  the  rear  HUMMV  wheel 

*  U3  =  the  road  input  for  the  howitzer  wheel. 


6.  Assembling  Lagrange's  Equation 

Combining  the  parts  of  Lagrange's  equations  that  were 
.developed  in  the  above  paragraphs  results  in  a  system  of  equations 
of  motion  in  the  form: 

(M]{x}  +  [C(x)){x}  +  [K]{x}  =  (f] 

where: 

[M]  =  the  7  X  7  system  mass  matrix 

[C(x)]  =  the  7  X  7  system  damping  matrix  that  is  a  function  of 
whether  the  shock  absorbers  are  in  tension  or  compression 
[K]  =  the  7  X  7  system  stiffness  matrix 

{x},  {x},  &  {x}  =  7  X  1  column  vectors  that  contain  accelerations, 
velocities  and  positions  of  all  coordinates,  respectively. 

[f]  =  the  7x1  forcing  vector  that  results  from  the  road  inputs. 

The  [M],  lC(x)],  and  [K]  matrices  for  this  system  are  contained  in 
Appendix  B.  This  completes  the  development  of  the  heave  and 

pitch  EOMs.  The  resulting  EOMs  were  numerically  integrated  to 

determine  the  time  history  for  motions  of  the  component  masses  of 
the  system. 

C.  HUMMV/HOWITZER  MODEL  NUMERICAL  INTEGRATION 
The  following  section  describes  the  procedures  that  were  used 
to  numerically  integrate  the  previously  developed  equations  of 
motion. 

1.  Transformation  to  1st  Order  Equations 

In  order  to  numerically  solve  the  equations  of  motion,  the 
formulated  system  of  second  order  ordinary  differential  equations 
must  be  transformed  into  a  system  of  first  order  differential 


equations.  This  equation  transformation  was  performed  using  matrix 
algebra  and  results  in  the  following  expression  for  the  system  of  first 
order  differential  equations; 


[xdot] 


r  (0)  (11  1 

'x' 

0  ■ 

Li-M  'K]  (-M-'C(j)lJ 

x_ 

+ 

M‘f_ 

where: 


[xdot]  =  a  14  X  1  column  vector  of  which  the  first  seven  elements 
are  the  velocities  of  all  the  generalized  coordinates  and  the  last 
seven  elements  are  the  accelerations  of  the  generalized  coordinates. 


a  14  xl  column  vector  of  which  the  first  seven  elements 


are  the  positions  in  terms  of  the  generalized  coordinates  and  the  last 
seven  elements  are  the  velocities  of  the  generalized  coordinates. 


[0]  =  a  7  X  7  zero  matrix 


[I]  s=  a  7  X  7  identity  matrix 

[-M‘‘C(x)]  &[~M‘‘k]  =  these  are  7x7  matrices  that  result  from 
multiplying  the  negative  inverse  of  the  system  mass  matrix,  [M], 
times  the  damping  matrix,  [C],  and  the  stiffness  matrix,  [K], 
respectively 


ro 

Im'T 


=  a  14  X  1  vector  which  the  first  seven  elements  are  zero  and 


the  last  seven  elements  are  the  matrix  product  of  multiplying  the  inverse 
of  the  mass  matrix,  [M],  times  the  forcing  vector,  [f]. 


For  the  integration  the  system  was  started  from  rest,  thus,  all 
seven  of  the  system's  displacement,  {x},  and  velocity,  {x},  state  variables 
had  a  value  of  zero  for  initial  conditions. 


2.  Integration  Difficulties 


The  equations  of  motion  were  numerically  integrated  using  the 
Matlab  ODE  4S  integration  routine.  This  is  a  variable  time  step  4th  order 
Runge-Kutta  integration  routine.  Integrating  these  equations  of  motion 
using  ODE  45  was  complicated  by  the  non-linear  damping  of  the  shock 
absorbers  and  the  required  shutoff  of  the  tire  springs  when  they  leave 
contact  with  the  road  surface. 

When  the  damping  values  of  the  shock  absorbers  would  abruptly 
change  from  compression  to  tension  and  vice  versa,  there  would  be  a 

corresponding  abrupt  change  in  the  shock  absorber  force.  I>etermining 

when  the  shock  absorber  damping  value  should  change  is  built  into  the 

integration  code  by  checking  for  each  time  step  if  each  of  the  shock 
absorber  forces  are  compressive  or  tensile.  An  "if"  statement  ensures  that 
the  correct  damping  value  is  used  for  the  condition  that  exists.  The  abrupt 

changes  in  force  in  the  shock  absorbers  slowed  down  the  tolerance  based 
ODE  45  routine.  A  similar  situation  occurred  when  the  tire  spring  forces 

would  abruptly  change  to  the  respective  static  weights  when  contact  with 

the  road  was  lost  for  a  particular  wheel. 

To  increase  the  speed  of  the  ODE  45  routine  the  integration  was 
divided  into  one  second  time  segments  so  the  solution  matrix  that  contains 

the  state  variables  remains  relatively  small.  The  displacements  and 
velocity  states  at  the  end  of  the  previous  time  segment  are  used  as  the 

initial  conditions  for  the  next  one  second  time  segment. 

The  resulting  time  histories  for  the  displacements  and  velocities 

allowed  post-processing  to  be  completed  determine  accelerations,  forces. 
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and  stress  in  key  system  components.  The  computer  code  developed  for 
this  numerical  integration  is  contained  in  Appendix  C. 


3  5 


III.  ROADWAY  INPUTS 


This  chapter  explains  the  road  inputs  that  were  used  to  drive 
the  computer  model.  The  roadway  inputs  are  defined  in  terms  of  the 
vertical  change  upward  or  downward  of  the  values  of  the  road  inputs 
Ui,  U2,  and  U3  that  are  shown  on  Figure  2-10.  These  are  independent 
vertical  displacement  inputs  to  tires  for  each  of  the  wheels  in  the 
computer  model.  As  defined  in  Chapter  II,  there  is  one  input  for  the 
front  HUMMV  wheels,  (Ui),  a  second  input  for  the  rear  HUMMV 
wheels,  (U2),  and  a  third  input  for  the  howitzer  wheels,  (U3).  The 
model  can  be  driven  over  any  terrain  that  can  be  defined  in  terms  of 
the  road  inputs  Ui.  U2.  and  U3 .  The  road  inputs  are  functions  of  time 
which  allows  the  user  to  control  the  speed  at  which  the  model  is 
driven  by  using  the  time,  velocity,  and  distance  relationships.  For 
this  design  study  the  road  inputs  used  were  a  washboard  course  that 
was  like  one  of  the  courses  used  in  the  actual  test  at  Aberdeen 
Proving  Grounds,  MD  (see  Chapter  VI),  and  a  predicted  "worst  case" 
pothole  road. 

A.  WASHBOARD  COURSE 

The  washboard  course  that  was  used  to  drive  the  model  was 
constructed  to  be  similar  to  the  Six  Inch  Washboard  Course  at  APG. 
During  the  instrumented  test  at  APG  there  were  several  data  runs 
taken  on  the  Six  Inch  Washboard  Course  (again,  see  Chapter  VI). 
Comparing  the  model  results  for  this  type  of  washboard  course  to  the 
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actual  instrumented  data  from  the  APG  test  for  a  similar  washboard 
provided  a  means  of  validating  the  output  results  of  the  computer 
model. 

To  model  the  APG  Six  Inch  Washboard  Course  a  sinusoidal  input 
was  used  with  an  amplitude  of  three  inches  and  a  distance  from  peak 
to  peak  of  72  inches.  However,  because  the  distance  between  the 

front  HUMMV  wheels  and  both  the  rear  HUMMV  wheels  and  also  the 
howitzer  wheels  does  not  correspond  to  even  multiples  of  the  72  inch 
peak  to  peak  distance,  a  lag  had  to  be  built  into  the  input  sinusoid  for 
both  the  rear  HUMMV  wheels  and  the  howitzer  wheels.  This  lag  was 
necessary,  to  get  the  correct  timing  of  the  bump  wheel  inputs  for  all 
the  wheels  in  the  system. 

Additionally,  in  order  to  prevent  difficulties  in  the  integration 
routine  that  would  occur  if  the  wheel  inputs  were  all  turned  on  at 
the  start  of  the  integration  run,  the  rear  HUMMV  and  howitzer  wheel 
inputs  were  not  turned  on  initially.  These  inputs  were  turned  on  at 
the  time  when  the  respective  wheels  would  start  hitting  the 

washboard  course  bumps  if  the  HUMMV/howitzer  system  entered 
the  course  from  level  terrain.  In  other  words,  the  front  HUMMV 

wheels  entered  the  washboard  and  started  taking  road  input,  Ui, 
while  the  rear  HUMMV  wheels  and  the  howitzer  wheels  remained  on 
level  terrain  which  provided  zero  input  to  them.  Then  the  rear 
HUMMV  wheels  entered  the  washboard  and  started  taking  road 

input,  U2,  while  the  front  HUMMV  wheels  continued  to  have  the 
input,  Ui,  and  the  howitzer  wheels  continued  to  have  zero  input. 


Finally,  the  howitzer  wheels  entered  the  washboard  and  at  this  point 

all  three  wheels  had  their  separate  inputs  Ui,U2,  and  U3, 

This  sinusoidal  input,  for  which  the  wavelength  was  defined  in 
terms  of  horizontal  distance  between  peaks,  was  converted  into  a 
function  of  time  by  using  time,  velocity,  and  distance  relationships. 
To  calculate  the  time  between  successive  peaks  for  a  specified 
velocity  the  inputs  Ui,  U2,  and  U3  were  then  defined  in  terms  of  time 
which  allowed  the  computer  model  to  be  "driven"  at  any  desired 

speed.  Making  the  input  a  function  of  time  was  required  for  using 

the  time  based  ODE  4S  integration  routine  of  Matlab  to  solve  the 

equations  of  motion.  Figure  3-1  depicts  the  wheel  inputs  Ui,  U2,  and 
U3  for  the  washboard  course  being  driven  at  S.O  MPH. 


Figure  3-1  Six  Inch  Washboard  Course  Wheel  Inputs 


B.  "WORST  CASE"  POTHOLE  COURSE 

The  model  was  also  "driven"  over  a  trench  type  pothole  to 

simulate  some  of  the  most  severe  road  conditions  that  the  howitzer 
would  be  subjected  to  during  tactical  operations.  The  model  results 
for  A-Frame  force  from  this  terrain  were  used  to  predict  if  the  A- 

Frame  and  T-Bar  were  strong  enough  to  handle  this  "worst  case" 

terrain.  The  pothole  course  road  inputs,  Ui,  U2,  and  U3,  are 

determined  by  a  separate  computer  program  (see  Appendix  D)  that 
computes  the  changes  in  axle  height  as  each  of  the  HUMMV  wheels 
and  the  howitzer  wheels  roll  through  the  pothole.  The  road  inputs 

are  based  on  the  path  that  would  be  traced  out  by  the  center  of  a  17 
inch  radius  wheel  rolling  through  the  specified  pothole.  The 

computer  program  that  was  developed  to  determine  the  road  inputs 

allows  the  user  to  specify  the  desired  dimension  of  the  pothole  for 
driving  the  computer  model.  For  this  design  study  the  "worst  case" 
pothole  was  determined  to  be  20  inches  deep  and  30  inches  across 
with  85  degree  slopes  on  each  edge.  This  size  was  selected  because  it 
was  judged  to  be  the  largest  pothole  that  the  HUMMV  driver  would 
traverse  without  stopping  the  HUMMV.  Thus,  the  "worst  case" 

pothole  dimensions  were  selected  based  on  what  was  judged  to  be 
the  most  severe  terrain  the  HUMMV/howitzer  would  encounter 
while  traveling  during  tactical  operations.  The  "worst  case"  pothole 
selected  was  slightly  larger  than  the  dimensions  of  the  pothole  that 
was  used  in  an  informal  test  conducted  at  RTA  in  June  1992.  Figure 
3-2  shows  the  pothole  profile,  and  the  axle  path  traced  for  the 
howitzer  or  HUMMV  wheel  rolling  through  this  "worst  case"  situation. 
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IV.  MODEL  RESULTS 


This  chapter  contains  the  results  of  driving  the  computer  model 
over  the  terrain  that  was  explained  in  Chapter  III.  The  results  from 
driving  the  model  over  the  washboard  course  show  close  agreement 

actual  results  obtained  from  the  APG  test.  A  detailed  comparison  of 
the  model  results  to  the  actual  results  is  contained  in  Chapter  VI. 
The  model's  results  from  the  "worst  case"  pothole  were  used  to 

predict  the  maximum  value  of  A-Frame  force  for  the  howitzer.  From 

this  predicted  maximum  force  the  stress  in  the  A-Frame  and  T-Bar 
was  calculated  and  compared  to  the  yield  stress  for  the  material  (see 
Chapter  V). 

The  predicted  A-Frame  force  is  calculated  by  multiplying  the 
combined  spring  stiffness  of  the  A-Frame/T-bar  times  the  change  in 
length  of  A-Frame.  The  change  in  length  of  the  A-Frame  is  a 
function  of  the  difference  in  the  angular  displacements  of  the  gun 
tube  mass  and  the  howitzer  trail/saddle  mass  at  a  given  time.  The 
angular  displacements  of  the  gun  tube  and  the  trail  were  determined 
from  numerically  integrating  the  equations  of  motion. 

In  addition  to  determining  results  in  terms  of  A-Frame  force, 
the  system's  mass  and  stiffness  matrices  (see  Appendix  B)  were  used 
to  determine  natural  frequencies  and  mode  shapes  of  the  system. 
These  natural  frequencies  and  mode  shapes  were  used  to  gain 

physical  insight  into  the  system's  performance.  The  natural 
frequencies  and  mode  shapes  for  the  system  are  shown  in  Appendix 
E 
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The  Matlab  "FFT"  function  was  used  to  perform  the  frequency 
analysis  on  the  A-Frame  force  results  from  the  washboard  course  to 
determine  the  frequencies  at  which  the  A-Frame  force  occurred  in 
the  model  results.  This  frequency  analysis  a  used  1250  Hz  data 

sampling  rate.  This  sampling  rate  yielded  a  time  increment,  dt,  that 
equals  1/1 250  or  the  inverse  of  the  sampling  rate.  In  order  to  have 
a  frequency  increment  that  was  sufficiently  small,  8192  data  points 

were  used  for  the  "FFT".  The  record  length,  T,  equals  dt  times  the 

number  of  points  or  8192/1250.  The  frequency  increment,  df,  was 

the  inverse  of  the  record  length  or  1250/8192. 

A.  WASHBOARD  COURSE  RESULTS 

The  time  histories  of  the  computer  model's  predictions  for  A- 
Frame  force  from  the  Six  Inch  Washboard  Course  are  periodic  with  a 
period  that  corresponds  to  the  period  of  the  washboard  input.  But 
unlike  the  sinusoidal  washboard  input  the  A-Frame  force  output 
predictions  are  not  sinusoidal.  These  non-sinusoidal  A-Frame  force 

outputs  most  likely  results  from  the  effects  of  the  non-linear  shock 

absorber  damping  that  was  incorporated  into  the  model  (see  Chapter 

H). 

The  frequency  analysis  of  the  computer  model's  predictions  for 
A-Frame  force  shows  that  as  expected  the  largest  component  of 
output  force  occurs  at  the  same  frequency  as  the  washboard  input 

frequency.  In  all  cases,  there  is  also  a  significant  but  lesser  force 
component  that  occurs  at  a  frequency  that  is  twice  the  input 
frequency. 

The  computer  model  was  driven  over  the  six  inch  washboard 

course  at  speeds  of  3.75,  5.0  MPH,  and  5.25  MPH.  These  speeds 
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correspond  to  the  speeds  for  which  actual  data  was  taken  during  the 

APG  test.  The  A-Frame  force  levels  were  relatively  low  at  3.75  MPH. 

The  force  levels  in  the  A-Frame  reached  a  maximum  at  5.25  MPH 

because  this  speed  for  the  washboard  course  corresponds  to  a 
natural  frequency  of  the  system.  Driving  the  model  at  speeds  higher 
than  5.25  MPH  resulted  in  decreasing  A-Frame  force  levels.  The 
time  histories  and  frequency  analysis  of  the  model's  predicted 
results  for  A-Frame  force  are  shown  in  Figures  4-1  through  4-6. 

1.  Results  at  3.75  MPH 

The  time  domain  model  results  at  3.75  MPH  are  shown  in 
Figure  4-1.  The  predicted  A-Frame  forces  are  approximately  1500 
lbs  in  tension  and  1200  lbs  in  compression.  Figure  4-2  shows  the 
frequency  analysis  at  this  speed.  As  expected,  there  is  a  large  force 
component  at  the  input  of  3.75  MPH  or  a  frequency  of  0.92  Hz  and 
there  is  smaLcr  force  spike  at  twice  this  frequency,  1.84  Hz. 

2.  Results  at  5.0  MPH 

The  time  domain  results  are  shown  in  Figure  4-3.  The  A- 

Frame  force  levels  reach  2800  lbs  in  tension  and  4800  lbs  in 

compression.  Figure  4-4  shows  the  frequency  analysis  results.  Again, 
a  large  force  component  at  the  input,  5.0  Hz,  and  another  force  spike 

at  twice  this  frequency,  2.42  Hz. 

3.  Results  at  5.25  MPH 

The  time  domain  results  are  shown  in  Figure  4-5.  The  A- 

Frame  force  levels  reach  2900  lbs  in  tension  and  4700  lbs  in 

compression.  Figure  4-6  shows  the  frequency  analysis  results.  A 
large  force  component  exists  at  the  input  frequency,  1.28  Hz,  and 

smaller  force  component  exists  at  twice  this  frequency,  2.56  Hz. 
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Figure  4-1  Model  Results  for  A-Frame  Force  at  3.75  MPH  (Time 
History) 
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Figure  4-2  Model  Results  for  A-Frame  Force  at  3.75  MPH  (FFT) 


Figure  4-3  Model  Results  for  A-Frame  Force  at  5.0  MPH  (Time 


History) 
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Figure  4-4  Model  Results  for  A-Frame  Force  at  5.0  MPH  (FFT) 
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B.  NATURAL  FREQUENCIES  AND  MODE  SHAPES 

The  natural  frequencies  and  mode  shapes  of  the  system  were 
determined  by  using  the  mass  and  stiffness  matrices  that  were 

developed  in  Chapter  II.  The  eigenvalues  and  eigenvectors  of  the 
system  are  determined  from  the  matrix  that  results  from  multiplying 
the  inverse  of  the  mass  matrix,  [M*  >],  by  the  stiffness  matrix,  [K]. 

The  system's  natural  frequencies  are  the  square  root  of  the 
eigenvalues.  The  system's  mode  shapes  are  the  eigenvectors  that 
correspond  to  these  natural  frequencies.  Analysis  of  these  natural 
frequencies  and  mode  shapes  shows  that  the  howitzer's  bounce  mode 
occurs  at  a  natural  frequency  of  1.28  Hertz  which  corresponds  to  a 

speed  of  5.25  MPH  for  the  howitzer  traveling  over  the  washboard. 

By  providing  an  excitation  at  a  natural  frequency  of  the  system  (1.28 
Hz)  the  maximum  A-Frame  force  that  could  possibly  be  developed 
for  this  particular  washboard  input  resulted. 

C.  "WORST  CASE"  POTHOLE  COURSE  RESULTS 

The  model  results  from  the  "worst  case"  pothole  that  was 

identified  in  Chapter  III  show  that  maximum  A-Frame  force  occurs 

lbs  at  15  MPH.  At  speeds  slower  or  faster  than  this  speed  the  model 

results  showed  lower  levels  of  A-Frame  force  for  this  pothole.  Figure 

4-7  shows  the  A-Frame  force  for  the  "worst  case"  pothole  at  15  MPH. 

There  are  spikes  in  the  A-Frame  force  when  the  wheels  of  the 

system  hit  the  pothole.  On  Figure  4-7,  at  0.1  seconds  the  front 

HUMMV  wheels  hit,  at  0.6  seconds  the  rear  HUMMV  wheels  hit,  and 

at  1.2  seconds  the  howitzer  wheels  hit  the  pothole.  When  the 
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howitzer  wheels  hit  at  IS  MPH  the  maximum  A-Frame  force  is  8,500 
lbs. 


A-Frame  Force 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 


Time  (sec) 

Figure  4-7  Model  Results  of  A-Frame  Force  for  "Worst  Case"  Pothole 
at  15  MPH 


V.  STRUCTURAL  ANALYSIS 


This  chapter  explains  the  structural  analysis  that  was 
performed  on  the  howitzer  in  order  to  determine  both  the  force 
levels  in  the  A-Frame  and  the  stress  on  the  howitzer's  components. 
To  determine  the  actual  force  levels  in  the  A-Frame,  the  A-Frame 
was  instrumented  with  strain  gages  as  shown  in  Figure  5-1.  This 
instrumented  A-Frame  was  then  calibrated  to  determine  a 
calibration  factor  which  would  correlate  strain  gage  readings  to  A- 
Frame  force  levels.  This  calibrated  A-Frame  was  used  during  the 
dynamic  test  (see  Chapter  VI)  at  Aberdeen  Proving  Grounds  (APG)  to 
measure  the  actual  force  levels  in  the  A-Frame. 

A.  A-FRAME  CALIBRATION 

The  primary  purpose  of  calibrating  the  A-Frame  was  to  develop 
a  reliable  means  to  correlate  dynamic  strain  gage  measurements  on 
the  legs  to  the  amount  of  force  that  causes  the  strain.  Additionally, 
the  strain  measurements  obtained  during  calibration  were  used  to 
determine  the  value  of  the  equivalent  linear  spring  stiffness  of  the 
A-Frame  for  the  computer  model.  During  the  calibration  extra  strain 
gages  were  attached  to  help  identify  levels  of  stress  at  key  locations 
other  than  the  locations  needed  for  the  force  versus  strain 
calibration.  These  extra  gages  provided  information  that  was  used  to 
form  conclusions  about  the  stress  in  the  A-Frame,  discussed  later. 


Figure  S-1  Location  of  A-Frame  Strain  Gages 


The  A-Frame  was  loaded  in  both  a  vertical  and  lateral  load 

configuration  to  determine  if  there  were  any  significantly  different 

effects  between  the  two  loading  configurations.  Figures  S-2  &  5-3 

show  the  test  setup  for  both  these  configurations. 

1.  A'Frame  Instrumentation 

To  conduct  the  calibration,  the  A-Frame  was  instrumented 

with  13  strain  gages  (Figure  S-1).  The  data  taken  from  the  pair  of 

single  strain  gages  located  on  the  bottom  of  each  leg  was  used  for  the 

development  of  the  calibration  curve.  The  pair  of  single  strain  gages 

on  the  top  of  each  leg  were  used  to  verify  that  loads  in  the  legs  were 
mainly  axial  loads  as  opposed  to  bending  loads.  The  single  gage  on 
the  A-frame  cross  support  brace  and  the  single  gage  on  the  attaching 
clamp  for  this  brace  were  used  to  identify  stress  levels  in  these 
members.  The  latter  two  strain  gages  were  used  because  there  were 
some  known  cases  where  the  cross  support  brace  clamps  had 

cracked,  and  it  was  unknown  whether  these  cracks  were  the  result  of 
high  load  levels  or  caused  by  manufacturing  defects. 

2.  Support  Equipment 

The  calibration  was  conducted  using  the  lab  facilities  of  the 
Department  of  Mechanical  Engineering  at  the  Naval  Postgraduate 

School.  To  hold  the  A-Frame  during  the  calibration  a  test  support 

structure  was  built.  The  main  objective  of  this  support  stand  was  to 

closely  match  the  mechanical  effects  that  occur  in  the  A-Frame  when 
it  is  loaded  while  connected  to  the  T-Bar.  The  test  support  used  was 
constructed  out  of  heavy  steel  and  used  actuti  eyebolts  from  a  Ml  19 
Howitzer  to  connect  the  bottom  of  the  A-Frame  legs  to  the  test 
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Figure  5-2  Vertical  Configuration  for  A-Frame  Test 


support  (Figure  S-4).  A  clamping  device  to  attach  the  apex  of  the  A- 
Frame  to  the  load  cell  was  also  designed  which  simulated  the  loading 

effects  of  the  actual  clamping  device  on  the  howitzer  (Figure  5-5). 

The  test  support  structure  achieved  the  goal  of  simulating  the  A- 
Frame  being  attached  to  the  T-Bar  of  the  howitzer. 

3.  Load  Cells 

The  vertical  loads  were  applied  using  a  Riehe  Testing 
machine  as  shown  in  Figure  5-2.  This  loading  device  was  used  to 
apply  vertical  loads  up  to  12,000  lbs.  The  lateral  loads  were  applied 
using  a  turn-buckle  that  was  connected  to  a  5,000  pound  rated 

capacity  dynamometer.  This  dynamometer  and  turn-buckle  system 
is  shown  in  Figure  5-3. 

4.  Calibration  Procedures 

Two  data  runs  were  recorded  for  both  the  vertical  and 

lateral  configurations.  A  data  run  consisted  of  loading  the  A-Frame 
in  increments  to  the  maximum  load  for  that  particular  data  run  and 
then  unloading  it  in  the  same  manner.  All  strain  gages  were  zeroed 
before  the  start  of  each  run.  At  each  load  increment,  strain  gage 

readings  were  recorded  for  all  of  the  strain  gages.  Tables  F-1 
through  F-4  contain  the  strain  gage  readings  and  load  data  obtained 
from  these  calibration  procedures. 
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5.  Calibration  Factor 

From  the  vertical  test  strain  gage  data  in  Appendix  F  the 
A-Frame  vertical  force  versus  strain  calibration  factor  was 
determined.  Since  the  A-Frame  vertical  load  test  was  a  static  test 
the  vertical  force  loading  at  the  apex  of  the  A-Frame  must  be 
balanced  by  the  sum  of  the  vertical  reaction  forces  in  the  legs  of  the 
A-Frame  (Figure  5-6). 


r 


where: 

F  =  the  applied  vertical  load  at  the  apex  of  the  A-frame 

^  horizontal  reaction  forces  where  the  A-frame 

connects  to  the  T-BAR 

Rvt  &  1^2  “  vertical  reaction  forces  where  the  A-frame 
connects  to  the  T-BAR _ 

Figure  5-6  Balance  of  Forces  on  the  A-Frame 
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The  reaction  forces  (Rv  i  &Rv  2)  at  the  bottom  of  each  leg  of 
the  A-Frame  were  determined  by  using  the  strain  gage  readings 
from  Gage  Pair  1  (gages  1  &  2)  and  Gage  Pair  4  (gages  9  &  10)  that 
are  identified  in  Figure  S-6.  The  average  value  of  the  strain  for 
either  of  the  strain  gage  pairs  was  multiplied  by  the  Modulus  of 
Elasticity  of  the  A-Frame  (10  x  10^  psi)  which  yields  the  axial  stress 
in  the  A-Frame  leg.  This  axial  stress  value  was  then  multiplied  by 
the  cross  sectional  area  of  the  A-Frame  legs  (0.6107  in^)  to 
determine  the  axial  force.  To  calculate  the  vertical  component  of  this 
axial  force  the  above  value  for  axial  force  was  multiplied  times  the 
cosine  of  29°.  A  sample  calculation  is  shown  below  using  the  10,000 
lb  load  increment  of  Vertical  Test  #1  (see  Table  F-1  for  strain  gage 
data): 

Strain  Gage  Pair  1  Avg  Strain  *  988  x  10"* 

Stress  =  Strain  x  Modulus  of  Elasticity 

-  (988  *  10'‘Xl0  *  10*psi)  -  9,880  psi 
»» Force  =  Stress*  Area  *  (9,880  psiXO.6107  in^) * 6,033  lbs 
.-.The  vertical  reaction  force:  R,,* (6,033  lbsXcos(28))  =  5,327  lbs 

Strain  Gage  Pair  4  Avg  Strain  “  1,095  *  lO"* 

Stress  "  Strain  x  Modulus  of  Elasticity 

=  (1,095  *10'^Xl0’‘10*psi)  =  1,095  psi 
=*  Force  *  Stress*  Area  « (10,950  psiXO.6107  in^) *  6,687  lbs 
.-.The  vertical  reaction  force:  R„2“ (6,687  lbsXcos( 28)) - 5,904  lbs 

T  Vertical  Forces  R„  +  =  5,327  +  5,904  =  11,231  lbs 


To  determine  the  calibration  factor  for  vertical  loading, 

the  percentage  difference  between  the  known  apex  load  and  the  sum 

of  the  calculated  vertical  force  components  at  the  bottom  of  the  A- 
Frame  was  determined.  The  average  percent  difference  for  all  the 
load  increments  is  equal  to  the  calibration  factor.  This  calibration 
factor  can  then  be  used  to  compare  measured  loads  in  the  A-Frame 
legs  to  the  applied  load  at  the  ^x  of  the  A-Frame.  An  average 
calibration  factor  of  0.88  was  determined  over  all  of  the  load 

increments  from  both  vertical  tests.  A  sample  calculation  of  the 

calibration  factor  for  the  10,000  lb  force  increment  of  Vertical  Test 

^1  (see  Table  F-1)  is  shown  below; 

Actual  Vertical  Force  lOtOOO  lbs 

Measured  Vertical  Force  "*U,23l  lbs 

CaUbration  Factor  -12:292-0.89 

11,231 

6.  A-Frame  Stiffness 

The  measured  strain  data  was  also  used  to  determine  a  n 

effective  linear  spring  stiffness  for  the  A-Frame.  To  determine  this 

stiffness  the  applied  vertical  load  for  a  speciHed  load  increment  was 
divided  by  the  corresponding  amount  of  vertical  deflection  in  the  A- 
Frame.  The  average  axial  deflection  in  the  A-Frame  legs  was 

calculated  from  the  average  strain  at  the  bottom  of  each  leg.  The 
average  axial  deflection  was  then  multiplied  by  the  cosine  of  28<*.  to 
calculate  the  amount  of  vertical  deflection.  A  sample  calculation  for 
the  10,000  lb  load  increment  of  Vertical  Test  #1  (see  Table  F-1)  is 
shown  below: 


60 


Strain  Gage  Pair  1  -  998;  Strain  Gage  Pair  4  »  1095 

QOfi  ^  109^ 

Ayg  axial  strain  * - - - -  1042  *  10"*;  Multii^y  by  leg  length 

Jm 

-•  Avg  axial  deflection  =  (28  in)(l042  *  10'*)  -  0.029  in 
Convert  to  vertical  deflection  =•  (0.029  inXcos(28))  =  0.0257  in 

.'.A-Frame  Stiffness  *  “  388,400  lbs/ in 

Q0257  in  _ 


7.  Observations 

From  analyzing  the  results  of  both  the  vertical  and  lateral 
A-Frame  loadings  the  following  conclusions  were  reached  regarding 
A-Frame  loading: 

*  The  strain  measurements  in  both  the  cross  member  brace  and 
the  attaching  clamp  for  this  brace  on  the  A-Frame  (strain  gages  3 
&  4,  in  Figure  5-1)  were  very  low  regardless  of  the  force  at  which 
the  A-Frame  was  loaded.  This  indicates  that  these  members  are 
not  critical  components  of  the  A-Frame  in  regards  to  failure. 

*  The  A-Frame  loads  are  carried  mostly  through  axial  loads  in  the 
A-Frame  legs  as  opposed  to  bending  loads.  This  means  that  the 
bottom  connection  joints  where  the  A-Frame  attaches  to  the  T-Bar 
act  as  pinned  joints  However,  because  there  was  some  bending  in 
the  A-frame  legs  at  lower  force  levels  the  bottom  joints  do  not  act 
as  perfect  pin  joints.  The  higher  force  levels  break  the  static 
friction  effects  in  the  joints  causing  the  joints  to  act  more  as 
perfect  pin  joints. 

B.  STRESS  ANALYSIS 

The  stress  analysis  focused  on  determining  the  highest  stress 
levels  in  the  T-Bar  and  the  A-Frame.  To  determine  if  mechanical 


failure  will  occur  due  to  exceeding  the  yield  stress,  the  stress  at  the 


probable  points  of  highest  stress  in  the  A-Frame  and  T-Bar  were 
compared  to  the  yield  stress  values  for  the  respective  materials. 

1.  T-Bar  Stress 

The  T-Bar  and  A-Frame  assembly  are  shown  below  in 

Figure  5-7,  Point  A  on  this  figure  is  identified  as  the  point  of  highest 
stress.  This  is  where  the  geometry  of  the  T-Bar  assembly  changes 
and  the  stronger  built  up  portion  of  the  T-Bar  ends.  To  simplify  the 

stress  analysis  the  T-Bar  is  modeled  as  a  cantilevered  beam.  Figure 
S-8  below  depicts  the  cantilevered  beam  modeling  of  the  T-Bar  as 

viewed  from  the  end  of  the  T-Bar  and  the  rear  of  the  howitzer.  The 
point  of  highest  stress  was  determined  to  be  at  the  comer  of  the  T- 
Bar  at  the  base  of  the  modeled  cantilevered  beam.  One  of  these 

comers  is  identified  as  Point  A  in  Figure  5-8. 


Figure  5-7  T-Bar  and  A-Frame  Assembly 
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Figure  5-8  T-Bar  Cantilevered  Beam  Model 

To  calculate  the  stress  at  Point  A,  the  normal  stress 

components  that  result  when  the  axial  force  in  the  A-Frame  leg,  P, 

was  resolved  into  components  in  the  x,  y,  and  z  directions  as  shown 

on  Figure  5-8.  The  stress  at  Point  A  is  the  sum  of  these  normal  stress 
components  using  the  principle  of  superposition  of  stress.  To  find 
the  normal  stress  components  first  the  A-Frame  leg  force,  P,  was 
projected  into  a  force  component  along  the  axis  of  the  T-Bar  in  the  y 
direction,  Py  =  sin(28)P,  and  a  second  force  component  in  the  xz 
plane,  Pxz=  cos(28)P.  The  Pxz  component  was  further  broken  up  into 
components  in  the  x  and  z  directions:  Px  =  cos(27)Pxz  and  Pz  = 

sin(27)Pxz-  To  calculate  the  normal  stress  that  was  caused  by  these 

force  components  Px,  Py.  and  Pz  well  known  strength  of  materials 

formulas  for  stress  in  a  beam  are  applied.  Therefore,  applying 

these  formulas  the  expression  for  the  maximum  stress  in  the  T-Bar 


where: 

Px,  Py,  and  Pz:  are  defined  in  the  above  paragraph. 

L=8  inches;  The  distance  from  the  base  of  the  cantilever  to  the 
point  of  application. 

A=.9240  in^:  The  cross  sectional  area  of  the  T-Bar. 

Ix=-8962  in^:  The  bending  moment  of  inertia  about  the  x  axis. 
Iz=.5979  in^:  The  bending  moment  of  inertia  about  the  z  axis. 

X  =1.300  in;  The  distance  in  the  x  direction  from  the  geometric 

center  of  the  T-Bar  to  the  outer  edge  of  the  T-Bar. 

z  =1.005  in:  The  distance  in  the  z  direction  from  the  geometric 

center  of  the  T-Bar  to  the  outer  edge  of  the  T-Bar. 


2.  A-Frame  Stress 

The  probable  point  of  highest  stress  in  the  A-Frame  will 
occur  around  the  0.37S  inch  diameter  transverse  circular  hole  that 
positions  the  cross  member  attaching  clamps.  This  hole  is  identified 
in  Figure  S-7  above.  At  the  location  of  this  hole  a  stress 
concentration  will  exist.  Figure  5-9  below  contains  the  empirical 
formula  from  page  601  of  [Ref.  4]  that  was  used  to  determine  the 
stress  concentration  in  a  pipe  section  with  a  transverse  hole. 

From  the  above  formula  the  stress  concentration  factor  for 
axial  loads  was  determined  to  be  3.73.  Additional  stress 
concentrations  due  to  the  effects  of  bending  and  torsional  loads  were 
neglected.  The  justification  for  neglecting  the  effects  of  bending 
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Figure  S'9  Stress  Concentration  Due  to  a  Circular  Hole  in  a  Pipe 


Section 


loads  was  that  both  the  A-Frame  calibration  results  (see  Appendix  F) 
and  the  results  of  the  instrumented  test  (see  Chapter  VI)  show  that 
the  joint  where  the  A-Frame  connects  to  the  T-Bar  acts  as  a  nearly 
perfect  pin  joint  at  higher  load  levels.  Any  stress  concentration 
reduction  due  to  the  spring  pin  that  is  driven  into  the  hole  were 
neglected  because  this  stress  reduction  was  assumed  to  be  small. 
Therefore,  maximum  total  stress  in  the  A-Frame  was  based  on  only 
the  axial  stress  component  in  the  legs  and  was  determined  using  the 
following  formula: 


o 


PK 

A 


(Eqn  5-2) 


where: 

P  =  the  axial  load  in  the  A-Frame  leg. 

K  =  3.73:  the  stress  concentration  factor  for  axial  loading  that  is 
caused  by  the  hole  in  the  A-Frame  leg. 

A  =  0.6107  in^:  the  cross  sectional  area  of  the  A-Frame  leg. 
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Without  this  hole  in  the  A-Frame  the  maximum  allowable 
stress  would  be  increased  by  the  amount  of  the  stress  concentration 

factor,  3.73.  This  hole  appears  to  only  aid  in  assembling  the  A-Frame 

and  serves  no  other  function.  Elimination  of  this  hole  would  make 
the  factor  of  safety  more  than  three  times  greater. 

3.  Worst  Case  Stress  and  Factors  of  Safety 

Using  Eqn  5-1  and  Eqn  5-2  above,  the  maximum  stress  was 
calculated  for  the  T-Bar  and  the  A-Frame.  The  maximum  predicted 
A-Frame  force  of  8,500  lbs  from  driving  the  computer  over  the 
"worst  case"  pothole  described  in  Chapter  II  was  divided  by  two  to 

calculate  the  for<.j  in  each  leg  of  the  A-Frame.  The  resulting  axial 
force,  P,  in  each  leg  of  the  A-Frame  equals  4,250  lbs.  From  Eqn  5-1 
the  resulting  stress  in  the  T-Bar  is  63,900  psi.  From  Eqn  5-2  the  A- 
Frame  stress  equals  26,000  psi. 

The  above  stress  levels  are  used  to  calculate  the  predicted 
factors  of  safety  for  the  T-Bar  and  the  A-Frame.  The  factor  of  safety 
is  calculated  by  dividing  the  yield  stress  of  the  component  by  the 
predicted  maximum  load  stress.  The  yield  stress  for  the  T-Bar  and 
A-Frame  are  116,000  psi  and  50,000  psi,  respectively.  Based  on  the 

maximum  predicted  stress  levels  given  above,  the  factor  of  safety  is 
1.81  for  the  T-Bar  and  1.92  for  the  A-Frame.  These  factors  of  safety 
show  that  both  the  T-Bar  and  the  A-Frame  are  strong  enough  to 
withstand  the  stress  from  the  worst  case  predicted  loads. 
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VI.  INSTRUMENTED  TEST 


An  instrumented  dynamic  test  was  conducted  on  the  Ml  19 
howitzer  at  Aberdeen  Proving  Grounds  (APG)  in  April  of  1993.  The 
howitzer  was  instrumented  as  described  below  in  paragraph  A.  The 
dynamic  testing  was  conducted  by  driving  the  howitzer  over  several 
of  the  APG  courses  at  various  speeds.  The  calibrated  A-Frame 
described  in  Chapter  V  was  used  to  measure  the  force  levels  in  the 
A-Frame  during  the  test.  In  most  respects,  the  test  results  for  A- 
Frame  force  are  similar  to  the  results  obtained  from  the  computer 
model.  However,  because  certain  physical  effects  were  not  included 
in  the  computer  model,  there  are  some  differences  between  the 
results  of  the  actual  test  and  the  results  from  the  computer  model. 
The  results  of  this  instrumented  test  confirmed  that  the  M 1  1 9 
Howitzer  could  be  safely  transported  in  the  firing  position  with  the 
A-Frame  attached. 

A.  INSTRUMENTATION 

For  this  test  the  howitzer  was  instrumented  with  20  strain 
gages,  two  tri-axial  accelerometers,  and  two  extensiometers.  This 
instrumentation  is  explained  in  the  following  paragraphs. 

1.  Accelerometers 

The  howitzer  was  instrumented  as  shown  in  Figure  6-1.  A 
tri-axial  accelerometer  was  placed  on  the  T-Bar  to  measure  the 
vertical,  longitudinal,  and  transverse  accelerations  of  the  howitzer 
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trail  mass.  This  accelerometer  was  not  placed  at  the  center  of 

gravity  of  the  trail  mass  because  of  difficulties  in  attaching  an 
accelerometer  at  this  point.  The  location  chosen  on  the  T-Bar 
provided  accelerations  on  the  same  order  of  magnitude  and 
representative  of  the  accelerations  at  the  center  of  gravity  of  the 
trail.  The  other  tri-axial  accelerometer  was  attached  on  the  elevating 
mass  near  its  center  of  gravity.  This  location  is  shown  in  Figure  6-2. 

2.  Single  Strain  Gages 

The  two  pair  of  single  strain  gages  on  the  bottom  of  each  leg 
(gage  pair  1  &  4  in  Figure  5-6)  of  the  previously  calibrated  A-Frame 
were  used  to  measure  the  strain  in  the  A-Frame  legs.  The  calibration 
factor  determined  in  Chapter  V  was  used  to  correlate  measured 
strains  in  the  A-Frame  legs  to  force  levels  in  the  A-Frame.  Single 
strain  gages  were  also  attached  on  the  top  and  the  bottom  of  the  T- 
Bar  on  both  sides  of  the  howitzer.  The  location  of  these  gages  is 
shown  in  Figure  6-3.  Because  the  highest  levels  of  stress  on  the  T- 
Bar  are  at  the  corner'^,  it  was  not  possible  to  place  the  gages  at  these 
points.  The  gages  were  placed  as  close  to  the  high  stress  points  as 
possible.  Even  though  these  gages  could  not  be  used  to  determine 
the  maximum  stress  on  the  T-Bar,  they  could  be  used  to  roughly 
judge  the  level  of  stress. 

3.  Other  Gages 

Strain  rosettes  were  placed  on  the  bottom  of  the  A-Frame 
legs  and  on  the  bottom  howitzer  trail  on  each  side  of  the  howitzer. 
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An  extensiometer  was  used  on  each  shock  absorber  to  determine  if 
the  shocks  were  extending  far  enough  so  that  the  stiff  rubber  bump 
stops  would  be  activated. 

B.  TEST  COURSES 

Several  of  the  standard  APG  test  courses  were  used  for  the 
instrumented  test.  These  courses  are  described  below. 

1.  Belgian  Block  Course 

This  course  is  designed  using  "Belgian  Blocks"  to  produce  a 
road  that  has  high  frequency  random  bumps.  The  magnitude  of  each 
bump  is  small.  This  course  is  shown  in  Figure  6-4.  During  the  test 
this  course  was  traversed  at  speeds  ranging  from  10  to  25  MPH. 


2.  Six  Inch  Washboard  Course 

This  course  has  six  inch  bumps  that  are  evenly  spaced  six 
feet  or  72  inches  apart  as  shown  in  Figure  6-S.  This  course  produces 
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significant  heave  and  pitch  effects  on  wheeled  vehicles  even  at 
relatively  slow  speeds.  During  the  test  this  course  was  traversed  at 
speeds  ranging  from  three  to  six  MPH. 


Figure  6-5  APG  Six  Inch  Washboard  Course 


3.  Cross  Country  Course  Number  1 

This  course  is  similar  to  a  relatively  smooth  but  winding 
trail  road.  This  course  was  traversed  at  approximately  25  MPH. 

4.  Cross  Country  Course  Number  3 

This  is  similar  to  a  rough  trail  road  which  has  random  large, 
medium,  and  small  sized  bumps.  This  course  was  traversed  at 
speeds  ranging  from  5  to  15  MPH. 
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C.  TEST  PROCEDURES 

Each  of  the  above  courses  were  traversed  by  the  howitzer  being 
towed  by  the  HUMMV  at  speeds  within  the  ranges  listed.  Data  was 
collected  on  all  the  channels  described  above,  using  a  sampling  rate 

of  2,500  Hertz.  Because  of  memory  limitations  in  the  data  collection 
equipment  only  approximately  30  seconds  worth  of  continuous  data 

could  be  stored.  For  each  course  several  data  runs  were  completed 

using  various  speeds.  For  instance,  on  the  Six  Inch  Washboard 
Course,  six  data  runs  were  taken  at  speeds  ranging  from  three  to  six 
MPH. 

The  speeds  for  each  data  run  on  the  Six  Inch  Washboard  Course 

and  Cross  Country  Course  Number  3  were  measured  using  the 
HUMMV's  speedometer,  because  these  courses  were  too  rough  to  use 

more  accurate  devices.  For  the  data  runs,  on  these  courses,  the 

driver  tried  to  hold  the  HUMMV  at  the  specified  speed  while  himself 
being  bounced  around.  Because  of  the  extremely  rough  road 

conditions  and  inaccuracies  in  the  speedometer,  the  actual  speeds 
ended  up  being  different  than  the  speeds  specified  in  the  test  plan. 
Consequently,  the  only  accurate  method  of  determining  the  actual 

speed  during  a  data  run  on  these  courses  is  from  the  video  tape 

taken  during  the  run. 

D.  TEST  RESULTS  AND  DISCUSSION 

The  Belgian  Block  and  Cross  Country  Course  Number  1  showed 
no  signiricant  levels  of  strain  or  acceleration  even  at  speeds  up  2  5 

MPH.  These  tests  simply  proved  that  the  howitzer  could  traverse 


these  courses  without  causing  damage  to  the  howitzer.  The  Six  Inch 
Washboard  Course  and  Cross  Country  Course  Number  3,  however,  did 
produce  significant  levels  of  strain  and  acceleration  on  the  howitzer 
components. 

The  test  results  for  the  Six  Inch  Washboard  Course  were 
analyzed  in  both  the  time  and  frequency  domains.  Data  runs  at 
speeds  of  3.75,  5.0  and  5.25  MPH  were  analyzed  for  the  Six  Inch 
Washboard  Course  and  compared  to  the  predicted  results  from  t  h  e 
computer  model.  The  most  severe  data  run  from  Cross  Country 
Course  Number  3  at  15  MPH  was  analyzed  in  the  time  domain.  The 

procedures  used  and  the  analysis  of  these  results  are  as  follows: 

1.  Time  History  Results 

The  time  history  for  vertical  "in  plane"  A-Frame  force  in 
each  case  was  produced  by  first  summing  the  "in  plane"  axial  forces 

from  both  of  the  A-Frame  legs.  The  same  calculation  procedures  as 
in  Chapter  V  were  used  for  determination  of  the  axial  forces  in  the 
legs  from  the  strain  gage  data  for  strain  gage  pairs  1  and  4  (see 
Figure  5-6).  The  vertical  force  was  determined  by  multiplying  the 
axial  force  times  the  cosine  of  28<^.  This  total  vertical  force  in  the  A- 
Frame  legs  was  multiplied  by  the  calibration  factor  of  0.88  that  was 
determined  in  Chapter  V.  The  A-Frame  lateral  forces  were 

calculated  by  determining  the  difference  in  the  horizontal  force 

components  of  the  axial  force  in  each  A-Frame  leg.  The  accelerations 

were  plotted  directly  from  the  measured  data. 

The  force  and  acceleration  time  histories  for  the  washboard 
course  (Figures  6-6  through  6-11)  show  that  high  frequency  spikes 


occur  on  the  washboard  course.  These  spikes  are  caused  the  braking 
system  being  activated  as  the  howitzer  went  over  the  bumps  for  all 

speeds.  Also,  for  the  S.2S  MPH  data  run  the  spade  hit  the  tops  of  the 
bumps.  These  effects  are  not  included  in  the  computer  model  and 

result  in  differences  between  the  model  results  and  the  actual 
results,  particularly,  for  the  S.2S  MPH  data  because  the  spade  was 
hitting. 

Time  histories  for  A-Frame  force  for  speeds  3.75,  5.0  &  5.25 
MPH  from  the  actual  test  (Figures  6-6,  6-8,  &  6-10)  are  similar  to  the 
time  histories  predicted  by  the  computer  model  for  the  same  speeds 

(Figures  4-1,  4-3  &  4-5).  The  high  frequency  force  spikes  do  not 

appear  in  the  model's  predicted  results,  however,  the  levels  of  the 

basic  force  response  in  the  actual  results  closely  compare  to  the  force 
response  of  the  model.  For  example,  the  actual  results  for  5.0  MPH 

(Figure  6-8)  has  spikes  that  reach  a  maximum  force  of  5,500  lbs.  But 
the  basic  force  response  on  average  reaches  approximately  2,400  lbs 

in  tension  and  4,300  lbs  in  compression.  The  corresponding  model 
prediction  (Figure  4-3)  shows  forces  of  2,800  lbs  in  tension  and 

4,800  lbs  in  compression.  A  similar  comparison  can  be  made  at  the 

3.75  and  5.25  MPH  speeds.  Although  the  actual  results  are  not 

exactly  the  same  as  the  model  predictions,  they  are  close  enough  to 

conclude  that  the  model  is  reasonably  representative  of  the  actual 
results. 

The  A-Frame  vertical  force  results  for  the  most  severe 

portion  of  the  Cross  Country  Course  Number  3  are  shown  in  Figure  6- 
12.  These  results  show  that  on  a  rough  trail  road  the  A-Frame  force 


Force  (lbs) 


is  not  as  great  as  the  force  that  was  caused  by  the  Six  Inch 
Washboard  Course  at  a  speed  of  5.25  MPH. 


Figure  6-6  A-Frame  Force  at  3.75  MPH  on  the  Washboard  Course 


Figure  6-7  Trail  Vertical  Acceleration  at  3.75  MPH  on  the 
Washboard  Course 
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Figure  6-8  A-Frame  Force  at  5.0  MPH  on  the  Washboard  Course 


79 


Accel  (Gs) 


Figure  6-9  Trail  Vertical  Acceleration  at  5.0  MPH  on  the  Washboard 


Course 
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Figure  6-10  A-Frame  Force  at 


1 


Force  (lbs) 


Figure  6-12  A-Frame  Force  at  15  MPH  on  Cross  Country  Course 
Number  3 
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2.  Frequencv  Domain 

The  frequency  domain  analysis  for  the  actual  A-Frame 
force  and  trail  acceleration  results  are  shown  in  Figures  6-13  through 
6-18.  The  same  procedures  using  the  "FFT"  function  in  Matlab  as 
were  used  in  Chapter  IV  for  the  model  results  were  used  to  perform 
the  frequency  analysis  on  the  actual  results.  The  time  increment 
used  was  1/1250  seconds  and  was  determined  by  using  every  other 
data  point  of  the  original  data  that  was  taken  at  2500  Hz.  8192 

points  were  used  to  achieve  an  adequate  frequency  discretization. 
The  record  length  was  8192/1250  and  frequency  increment  was 
1250/8192.  The  FFT  analysis  of  both  A-Frame  force  and  trail 
vertical  acceleration  for  each  data  run  clearly  shows  two  frequency 
spikes.  The  lowest  frequency  and  largest  spike  corresponds  to  the 
forcing  frequency  or  road  input  to  the  howitzer  tires.  The  second 
largest  spike  corresponds  to  twice  the  forcing  frequency. 

The  frequency  analysis  for  actual  A-Frame  force  for  speeds 
of  3.75,  5.0,  &  5.25  MPH  (Figures  6-13,  6-15,  &  6-17)  compares 

reasonably  well  to  the  corresponding  frequency  analysis  from  the 
model  results  (Figures  4-2,  4-4,  &  4-6).  For  example,  the  FFT  results 
from  the  actual  data  at  5.25  MPH  (Figure  6-17)  shows  a  large  force 

component  1180  lbs  at  the  road  input  frequency  (1.28  Hz)  and 
second  large  force  component  of  575  lbs  at  double  the  road  input 

(2.56  Hz).  Where  as  the  FFT  of  the  model  for  the  corresponding 
predicted  A-frame  force  (Figure  4-6)  shows  a  large  force  component 
of  1375  lbs  at  the  road  input  frequency  (1.28  Hz)  and  a  second  large 
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force  component  of  500  lbs  at  twice  the  input  frequency  (2.56  Hz). 
For  this  case,  the  total  amount  of  force  from  the  1.28  and  2.56  Hz 
spikes  for  the  actual  results  and  the  model  results  are  approximately 
the  same  (1765  lbs  versus  1875  lbs).  Although  FFT  results  for  the 


actual  data  and  the  model  predictions  are  not  exactly  the  same,  they 
are  close  enough  further  to  validate  the  model's  predicted  results. 


Vertical  Force  FFT  from  APG  Test  (3.75  MPH) 


Figure  6-13  FFT  of  A-Frame  Force  at  3.75  MPH  on  the  Washboard 
Course 


85 


Accel  (Gs) 


86 


Force  (lbs) 


Figure  6-15  FFT  of  A-Frame  Force  at  5.00  MPH  on  the  Washboard 


Course 
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Figure  6-17  FFT  of  A-Frame  Force  at  5.25  MPH  on  the  Washboard 
Course 
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Vertical  Trail  Accel  FFT  from  APG  Test  (5.25  MPH) 


Figure  6-18  FFT  of  Trail  Vertical  Acceleration  at  5.25  MPH  on  the 
Washboard  Course 
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3.  Other  Results 

The  strain  reading  from  the  strain  rosettes  that  were 

attached  to  the  bottom  the  howitzer's  trail  showed  low  strain 
readings  for  all  courses.  These  low  strain  readings  clearly  showed 
that  transporting  the  howitzer  in  the  firing  position  with  the  A- 
Frame  attached  produces  insignificant  levels  of  stress  in  the  howitzer 
trail. 

Lateral  forces  in  the  A-Frame  were  also  determined  from 
the  actual  test  data.  Figures  6-19  through  6-21  show  the  lateral 

force  results  for  3.75,  5.0  &  5.25  MPH,  respectively.  The  maximum 
lateral  forces  were  less  than  ten  percent  of  the  maximum  vertical 
forces  for  all  runs  on  the  washboard  course.  For  the  rough  terrain 
Cross  Country  Course  Number  Three  the  maximum  lateral  force 
(Figure  6-22)  was  twenty  percent  of  the  maximum  vertical  force. 
These  relatively  low  values  of  lateral  force  help  to  show  the  validity 

of  the  modeling  technique  used,  since  the  computer  model  only 
predicts  the  vertical  forces  in  the  A-Frame  and  assumes  that  the 

lateral  forces  would  be  comparatively  low. 
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Figure  6-20  Lateral  Force  from  the  APG  Washboard  Course  at  5.0 


MPH 
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Figure  6-22  Lateral  Force  from  Cross  Country  Course  Number  Three 
at  15  MPH 
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VII.  CONCLUSIONSIRECOMMENDATIONS 


A.  CONCLUSIONS 

The  computer  model  developed,  produced  output  results  that  are 
reasonably  representative  of  the  actual  results  for  the  washboard 
type  road  input  used  at  Aberdeen  Proving  Grounds.  Therefore,  the 
model  is  validated  and  can  be  used  to  predict  "worst  case"  loading  on 
the  howitzer's  critical  parts.  These  "worst  case"  predictions  for  a 
trench  type  pothole  show  the  Ml  19  Howitzer  can  be  safely  towed  in 
the  firing  position  with  the  A-Frame  attached. 

As  an  additional  result  of  this  study,  we  believe  that  the  T-Bar  is 
the  most  critical  component  on  the  howitzer  in  regards  to  mechanical 
failure  when  towing  in  the  firing  position.  The  predicted  factor  of 
safety  for  the  T-bar  based  on  the  "worst  case"  pothole  terrain  is  1.81. 
Due  to  the  high  cost  of  repair,  making  the  T-Bar  out  of  thicker  steel 
would  provide  an  increased  factor  of  safely. 

The  predicted  factor  of  safety  for  the  A-Frame,  from  the  same 
"worst  case”  terrain,  is  1.92.  This  indicates  that  the  A-Frame  is 
strong  enough  to  withstand  the  worst  case  expected  loads.  However, 
the  A-Frame  is  made  out  of  aluminum  and  aluminum  does  not  have 
a  fatigue  limit.  Thus,  eventually  the  A-Frame  is  expected  to  fatigue. 
The  A-Frame's  fatigue  life  could  be  greatly  extended  by  eliminating 
the  holes  in  its  legs. 
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The  strain  reading  taken  from  the  trail  rosettes  during  the 
instrumented  test  at  APG  show  that  the  stress  in  the  trail  is 
insignificant  during  towing  in  the  A-Frame  position. 

B.  RECOMMENDATIONS 

The  US  Army  should  remove  the  towing  restrictions  and  allow 
the  howitzer  to  be  towed  in  the  firing  configuration  with  the  A- 
Frame  attached  during  tactical  operations. 

The  T-Bar  should  be  made  out  of  thicker  steel  to  provide  added 
assurance  that  it  will  not  fail. 

The  A-Frame  should  be  re-configured  for  manufacture  without 
the  hole  in  its  legs. 
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APPENDIX  A 


SYSTEM  PARAMETER  VALUES 

The  following  table  defines  the  HUMMV/Howitzer  system 
parameters  using  an  inches  (in),  pounds  (lbs),  and  seconds  (sec) 
system  of  units.  Gravity  (G)  is  approximated  to  be  equal  to  386.4 
in/sec^ .  The  nomenclature  for  the  parameters  that  are  identified  in 
this  Appendix  is  the  same  nomenclature  that  is  used  throughout  this 
paper.  Refer  to  Figures  2-1  and  2-2  to  identify  the  parameters 
defined  below.  The  values  listed  in  the  table  incorporate  the  effects 
of  both  sides  for  those  components  that  are  identical  on  the  right  and 
left  side  of  the  HUMMV  or  Howitzer.  For  example,  the  spring 
stiffness  for  each  of  the  front  HUMMV  tires,  Ki,  is  1323.0  Ibs/in.  This 
value  is  multiplied  by  two  to  account  for  the  combined  effect  of  both 
front  HUMMV  front  tires.  The  value  listed  in  the  following  table  for 
Ki  is  2646  Ibs/in. 
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TABLE  A-1  COMPUTER  MODEL  PARAMETER  VALUES 


Symbol 

Value/units 

Ki 

2646.0  Ibs/in 

Kz 

1908.0  Ibs/in 

Js 

3544.0  Ibs/in 

K* 

5040.0  Ibs/in 

K, 

2000.0  Ibs/in 

K7 

1056.0  ibs/in 

mm 

200,000.0  ibs/in 

Bic»l71.6  IbS’sec/in 
Bu«308.4  Ibs-sec/ln 

Bz 

B2c«199.2  Ibs-sec/in 

Bs 

1  -fiS,  1 

B 

250,000  Ibs-in-sec/rad 

M, 

1.294  Ibm 

M2 

1.294  Ibm 

M3 

21.22  Ibm 

M4 

1.035  Ibm 

Ms 

9.834  Ibm 

Me 

5.952  Ibm 

h 

52,680.0  Ibs-in-sec^ 

How  Oetermined/Notes 


Reference  1/HUMMV  front  tires  stiffness 


Reference  1  /HUMMV  front  suspension  stiffness 


Reference  1/  HUMMV  rear  tires  stiffness 


Reference  1  /  HUMMV  rear  suspension  stiffness 


Estimated /Howitzer  tires  stiffness 


ererence  Z/Hodewd  Howitzer  vertical 
stiffness 


Calculated/Combine  A-Frame  &  T-Bar  stiffness 


Reference  1/Bic  is  for  compression  &  Bit 
for  tension.  Values  are  20%  more  than  Ref  1 


Reference  1/  Bzc  is  for  compression  &  Bzt 
for  tension.  Values  are  20%  more  than  Ref  1 


Calculated  from  RIA  data/  Bsc  is  for 
ion  &  Bst  for  tension. 


Estimated/  Torsional  Damping  in  elevating  clutch 


Reference  1/  Includes  mass  of  HUMMV  rear 
wheels,shocks  and  suspension  springs 


Refwencel/  HUMMV  Body 


Reference  2/Includes  howitzer  wheels, 


R  m*  !•  Wr-  \  i  iili  ir 


Reference  2/  Howitzer  trail  &  saddle  mass 


Reference  2/  Howitzer  elevating  mass 


Reference  1/  Moment  of  inertia  about  HUHMVCG 
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TABLE  A-1  COMPUTER  MODEL  PARAMETER  VALUES  (CONT) 
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APPENDIX  B 


MATRICES  FOR  EQUATIONS  OF  MOTION 
This  Appendix  contains  the  mass,  damping,  and  stiffness 
matrices  for  the  equations  of  motion  that  were  developed  in  Chapter 
II.  The  mass  system  matrix,  [M],  is  as  follows; 


[mi,  0  0  0  0  0  0  1 

lo  m,,  0  0  0  0  0  I 

I  I 

0  0  mjj  0  m,s  mjs  m,, 

jo  0  0  m^4  0  0  0  j 

jO  0  mj3  0  m„  m„  m,,  j 

0  0  m„  0  m„ 

[o  0  m73  0  m^j  m,,;  m77  J 

where: 

ill  1 1  *  III22  *  ^^2 

m,,  -  MjLj  +  MjL, 

®J7  *  ^6^CG  COS(0q) 
mjj  =  MjLj  +  MjLj 
mjj  =  M5L3LJ1  +  M(L3L,0 


m,«  - 12  +  M5L,,"  + 
®73  “  ^6^CO  COS(0q) 

*^76  *  MgL(.gLjgCOs(0o) 


m3j-M3  +  M5+M* 

mjj^MjLji  +  M^L,( 
m44  * 

mjs  *  Ii  +  MjL,^  +  MjL,^ 
®57  ~  ^6^CG^3 

m(5*MjL3Lii  +  M0L3L,6 
®67  “  ^6U{.QLjgCOS(00) 

®75  *  ^6^CG^3  Ct)S(0Q) 

®77  “  ^3  ^6^CG^ 


Note  due  to  symmetry 
equal: 

mj3  -  m35 

n>«  “  “56 


of  [M]  the  following  matrix 

***43  “  ***36  ***73  “  ***37 

***7S  "  ***57  ***76  *  ***67 


elements  are 
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The  damping  matrix  for  the  system  is  dependent  on  whether 
the  shock  absorber  forces  are  tensile  or  compressive  at  a  given  time. 
The  tensile  or  compressive  state  for  each  of  the  shock  absorbers  is 
determined  by  the  relative  velocities  of  the  end  connection  points  for 
the  shock  absorber.  The  damping  matrix,  [C(x)],  is  as  follows: 


fCn 

0 

C.3 

0 

Cl5 

0 

0  1 

lo 

^23 

0 

^25 

0 

0 

1 

^32 

C33 

C34 

C35 

^36 

0 

0 

0 

®43 

C44 

^45 

^46 

0 

1  ®51 

C52 

C53 

^54 

C55 

C56 

0 

1® 

0 

^63 

^64 

^65 

C66 

C67 

[o 

0 

0 

0 

0 

C76 

^7j 

where: 


u  “B, 

^13  “  “Bj 

Cj5  “  B,Lj 

22  “  B2 

C23  =  -B2 

Cjs  *  “B2L2 

31  "-B, 

C32  ”  “Bj 

C33  =  Bj  +  B2  +  B; 

34  “  “B3 

C35  BjL,  +  B2L2  +  B3L3 

36  “  B,L,j 

^43  “  ~B3 

C44  =  B3 

4S  “  “BjLj 

®46  “  ~BjLi2 

^51  “  BjLj 

52  “  “BjLj 

C53  =  “B,Lj  +  BjLj  +  B3L3 

54  “  “BjLj 

C55  *  B,Lj  +  BjLj  +  B3L3 

56  “  BjLjLjj 

^63  “  ^^12 

^64  ”  ”B3L,2 

65  “  ^^3^2 

^66  “  +  B4 

^67  “  “B4 

76  “  “B4 

C77  “  B4 

Note  due  to  symmetry  of  [C(x)]  the  following  elements  are  equal: 


31  “  ^13 

Cj2  "  C23 

^51  ^15 

*^43  *^34 

52  “  ^25 

C53  *  C35 

*'54  “  C45 

*^63  “  *^36 

64  “  ^4^ 

*'65  *  ^56 

^6  ”  *'67 

102 


The  system  stiffness  matrix,  [K],  is  as  follows: 


fk,, 

0 

k„ 

0 

k,s 

0 

0 

0 

^22 

k23 

0 

kjs 

0 

0 

kj. 

kj2 

kss 

k34 

k3S 

k36 

0 

0 

0 

k43 

k44 

k4S 

k46 

0 

kj, 

k52 

k53 

k54 

ks5 

ks6 

0 

0 

0 

k63 

k64 

kes 

k66 

k67 

Lo 

0 

0 

0 

0 

k76 

k?? 

ku“K.  +  K, 

ki3  “  ^2 

k,5=K2L, 

k22  -  K3  +  K, 

k  =  -K 

*■23  JV4 

k,5»-K4L2 

M 

1 

II 

k32  ”  -K4 

k„  =  K2  +  K,+K. 

k34  “  -K7 

k33  =  -K2L,  +  K4L2  + 

K7L3 

ks*  ”  K7L12 

k43  *  -K7 

k44  “  Kj  +  K7 

k^j  ™  ”1^71^3 

k46  “  “®^7^2 

k5,“K2L. 

kjj  “  ~K4L2 

k33  -  -K2L,  +  K4L2  + 

K7L3 

k54  “  “K7LJ 

k55  -  K2L,^  +  K4L2'  + 

K7L,^ 

ks6  “  K7L5Lj2 

k63  “  ^7^-12 

k64  “  ”^7^2 

k«5  “  ^7k'3Li2 

kj7  “  —  Ky^Ly^L|3COS(6Q) 

k76  “ 

k77  ”  Ky^Ly^Lj3COS(0Q) 

Note  due  to  symmetry 

of  [K]  the  following  matrix 

elements  are  equal: 

ksi  “  kj3  kj2  =  k23 

kji  =  k,5  k43  =  k34 

k52  *  kjs  k53  =  kjs 

ks4  *  k45  kjj  “  kjg 

k64  “  k4g  kg5  =  kjg 

k76  “  k57 
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APPENDIX  C 


NUMERICAL  INTEGRATION  COMPUTER  CODE 


This  Appendix  contains  the  computer  code  used  to  integrate 
the  equations  of  motion  using  the  ODE  45  integration  routine  in 
Matlab.. 

1.  THE  ODE  45  INTEGRATION  FUNCTION 

function  xdot=mateqn(t,x) 

%  Define  gravity 
G=386.4; 

%  Define  HUMMV  Parameters: 

Ll=83.8;  L2=46.7;  L3=75.7;  L7=130.5;  L8=294.5; 

Kl=2»1323,0;  K2=2*954.0;  K3=2*1777.0;  K4=2*2520.0; 

K5=1000000.0; 

B1C=2*1.2*71.5;  B1T=2*1.2*128.25;  B2C=2*  1.2*83.0; 

B2T=2*  1.2*200.0; 

Wl=2*250.0;  W2=2*250.0;  W3s8200.0; 

M1=W1/G;  M2=W2/G;  M3=W3/G; 

11=52680;  RW=5400;  FW=3800; 

%  Define  Gun  Parameters: 

L4=123.4;  L5=11.6;  LA=45.0;  LS=42.6;  LCG=32.5; 

LI  1=104.4;  L12=135.0;  L13=38.0;  L14=17.0;  L15=9.70; 

L16=110.3; 

K6=2*  1000.0;  K7=2*528;  KA=200000;  KS=300; 

B3C=2*20.0;  B3T=2*66.0;  B4=.5*250000.0; 

W4=2*200;  W5=1500;  W6=2300;  M4=W4/G;  M5=W5/G; 

M6=W6/G; 

12=10000;  13=15000; 

Q0=12.0*pi/180;  PHlA=51.0*pi/180;  PHIS=14.0*pi/180; 

AW=3800;  TB=88000;  %  Define  Mass  matrice 

%  Define  the  mass  matrix 
M=[zeros(7,7)]; 
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M(1,1)=M1;  M(2,2)=M2;  M(3,3)=M3+M5+M6;  M(3,5)=M5*L3+M6’‘L3; 

M(3,6)=M5*L1 1  +M6*L16;  M(3,7)=M6*LCG*cos(Q0); 

M(4,4)=M4;  M(5,3)=M5*L3+M6*L3;  M(5,5)=I1+M5*L3''2+M6''2; 

M(5,6)=M5*L3*L11+M6*L3*L16;  M(5,7)=M6*L3*LCG*cos(Q0); 

M(6,3)=M5*L11+M6*L16;  M(6,5)=M5*L3*L1 1+M6*L3*L16; 

M(6,6)=I2+M5*L1P2+M6*L16''2+M6*L14'^2; 
M(6,7)=M6*L16*LCG*cos(Q0); 

M(7,3)=M6*LCG*cos(Q0);  M(7,5)=M6*L3*LCG*cos(Q0); 

M(7,6)=M6*L16*LCG*cos(Q0); 

M(7,7)=I3+M6*LCG^2; 


%  Define  Stiffness  Matrix 
K=[zeros(7,7)l; 

K(1,1)=K1+K2;  K(l,3)=-K2;  K(1,5)=K2*L1; 

K(2,2)=K3+K4;  K(2,3)=-K4;  K(2,5)=-K4*L2; 

K(3,l)=-K2;  K(3,2)=-K4;  K(3,3)=K2+K4+K7;  K(3,4)=-K7;  K(3,5)=- 
K2*L1+K4*L2+K7*L3;  K(3,6)=K7*LI2; 

K(4,3)=-K7;  K(4,4)=K6+K7;  K(4,5)=-K7*L3;  K(4,6)=-K7*L12; 

K(5,1)=K2*L1;  K(5,2)=-K4*L2;  K(5,3)=-K2*L1+K4*L2+K7*L3;  K(5,4)=- 

K7*L3; 

K(5,5)=K2*L1''2+K4*L2''2+K7*L3''2;  K(5,6)=K7*L3*L12; 

K(6,3)=K7*L12;  K(6,4)=-K7*L12;  K(6,5)=K7*L3*L12; 

K(6,6)=K7*L12''2+KA*LA*L13*cos(Q0); 

K(6,7)=-KA*L13*LA*cos(Q0); 

K(7,6)=-KA*L13*LA*cos(Q0); 

K(7,7)=+KA*L13*LA*cos(Q0); 

%  Define  and  velocities: 

VHF=x(10)  -  (Ll*x(12));  VHB=x(10)  +  (L2*x(12)); 

VGB=x(10)  +  (L3*x(12))  +  (L12*x(12)); 

%  Define  Road  Inputs: 

VMPH=5.0;  VEL=VMPH*17.6;  PD=72A^EL; 

Ul=3*(sin(2*pi*t/PD)); 

if  t  <  L7A^EL,  U2=0;else,  T2=t-(L7A^EL); 

U2=3*(sin(2*pi*T2/PD));  end 

if  t  <  LSA^EL,  U3=0;else,  T3=t-(L8A^EL); 

U3=3*(sin(2*pi*T3/PD));  end 
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%  Determioe  the  shock  absorber  damping  constants  based  on 
compression  or  tension 


%  HUMMV  Front  shock 

if  (VHF  -  x(8))  >  0.0,  B1=B1T;  else,  B1=B1C;  end 
%  HUMMV  Rear  shock 

if  (VHB  -  x(9))  >  0.0,  B2=B2T;  else,  B2=B2C;  end 
%  GUN  shock 

if  (VGB  -  x(ll))  >  0.0,  B3=B3T;  else,  B3=B3C;  end 

%  Define  Damping  Matrix 
B=[zeros(7,7)]; 

B(1,1)=B1;  B(1,3)=-B1;  B(1,5)=B1*L1; 

B(2,2)=B2;  B(2,3)=-B2;  B(2,5)=-B2*L2; 

B(3,1)=-B1;  B(3,2)=-B2;  B(3,3)=B1+B2+B3;  B(3,4)=-B3;  B(3,5)=- 

B1*L1+B2*L2+B3*L3;  B(3,6)=B3*L12; 

B(4.3)=-B3;  B(4,4)=B3;  B(4,5)=-B3*L3;  B(4,6)=-B3*L12; 

B(5,1)=B1*L1;  B(5,2)=-B2*L2;  B(5,3)=-B1*L1+B2*L2+B3*L3; 

B(5,4)=-B3*L3;  B(5,5)=B1*L1''2+B2*L2''2+B3*L3''2; 

B(5,6)=B3*L3*L12; 

B(6,3)=B3*L12;  B(6,4)=-B3*L12;  B{6,5)=B3*L3*L12; 

B(6,6)=B3*L12''2+B4; 

B(6,7)=-B4;  B(7,6)=-B4;  B(7,7)=B4; 

%  Define  Forcing  Vector 
F=[zeros(7,l)]; 

F(1,1)=K1*U1;  F(2,1)=K3*U2;  F(4,1)=K6*U3; 

%  CHECK  IF  FRONT  HUMMV  WHEEL  LOSES  CONTACT 
if  K1*(U1  -  x(l))  <  -FW,  F(2,1)=-FW;  K(1,1)=K2;  end 

%  CHECK  IF  REAR  HUMMV  WHEEL  LOSES  CONTACT 
if  K3*(U2  -  x(2))  <  -RW,  F(2,1)=-RW;  K(2,2)=K4;  end 

%  CHECK  IF  HOWITZER  WHEEL  LOSES  CONTACT 
if  K6*(U3  -  x(4))  <  -AW,  F(4,1)=-AW;  K(4,4)=K7;  end 

%  transform  into  system  of  1st  order  ode's 

id=[l  zeros(l,6);0  1  zeros(l,5);0  0  1  zeros(l,4);zeros(l,3)  1  zeros(l,3)]; 
id=[id;zeros(l,4)  1  0  0;zeros(l,5)  1  0;zcros(l,6)  1]; 
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a=l2eros(7,7)  id;  -inv(M)*K  -inv(M)*Bl; 
fl=inv(M)*  F; 
f=[[zeros(7,l)];  fl]; 

%integrate 
xdot=a*x  +  f; 


2.  THE  DRIVER  CODE  FOR  ODE  45 

%  Set  the  Initial  Condition 
x01=[zeros(14,l)]; 

%  Integrate  in  one  second  increments 
to  1=0.00;  tf 1=0.999;  %x01=xn(length(xn),:); 
[tl,xl]=ode45('mateqn',t01,tfl,x01,le-3); 
t02=1.0001;  tf2=1.999;  x02=xl(length(xl),:) 
[t2,x2]=ode45('mateqn',t02,tf2,x02,le-3); 
t03=2.0001;  tf3=2.999;  x03=x2(length(x2),:) 
[t3,x31=ode45('mateqn’,t03,tf3,x03,le-3); 
t04=3.0001;  tf4=3.999;  x04=x3(length(x3),:) 
[t4,x4]=ode45('mateqn',t04,tf4,x04,le-3); 
t05=4.0001;  tf5=4.999;  x05=x4(length(x4),:) 
[t5,x5]*ode45('mateqn',t05,tf5,x05,le-3); 
t06=5.0001;  tf6=5.999;  x06=x5(length(x5),:) 
[t6,x6]=ode45('mateqn',t06,tf6,x06,le-3); 
%t07=6.0;  tf7=6.999;  x07=x6(length(x6),:); 
%[t7,x7]=ode45('mateqn',t07,tf7,x07,le-3); 
%t08=7.0;  tf8=7.999;  x08=x7(length(x7),:); 
%[t8,x8]=ode45(’mateqn',t08,tf8,x08,le-3); 

%  Combine  the  solution  state  matrix 
xn=[x  1  ;x2;x3  ;x4  ;x5  ;x6] ;  %  x7;x81 ; 
tn=[tl;t2;t3;t4;t5;t6];%t7;t8]; 
save  ghl  xn  tn 
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APPENDIX  D 


POTHOLE  COURSE  COMPUTER  CODE 

This  Appendix  contains  the  computer  code  that  was  developed 
to  model  the  terrain  inputs  for  a  road  with  a  trench  type  pothole.  This 
program  calculates  and  plots  a  trace  of  a  pothole  proHle  and  the  axle 
path  followed  by  circular  wheel  rolling  through  this  type  pothole.  This 
program  works  for  any  pothole  that  can  be  model  with  a  flat  bottom, 
and  sides  that  can  be  approximated  as  straight  lines.  The  portion  of 
the  code  shown  below  only  works  for  potholes  that  have  an  entrance 
height  that  is  equal  to  the  exit  height.  Additional  code  to  determine 
the  terrain  inputs  for  a  any  type  of  pothole  is  on  file  at  the  Naval 
Postgraduate  School.  Computer  code  was  also  developed  to  model 
bump  inputs,  such  a  log  in  a  road  and  maintained  on  file  at  the  Naval 
Postgraduate  School. 

COMPUTER  CODE: 

function  [x,y,ya,u]=pothole(xent,depth,thetal,xbot,theta2,farht,xtrail); 

%  The  call  statement  for  this  program/function  is  as  %  follows: 

%  [x,y,ya,u]=pothole(xent,depth,thetal  ,xbot,theta2,farht,xtrail) 

%  Input  parameters  used  in  the  call  are  defined  as  follows: 

%  depth  -  the  vertical  distance  from  the  bottom  to  the  level 
%  entrance  surface. 

%  theta  1  -  the  angle  of  the  decline  on  the  entrance  side,  measured 

%  from  an  extension  of  the  level  entrance  surface  to  the  decline. 

%  xbot  -  the  horizontal  distance  across  the  bottom. 

%  theta2  -  the  angle  decline  on  the  exit  side,  measured  like  theta  1. 

%  farht  -  the  vertical  distance  from  the  bottom  to  the  level  exit 

%  surface. 

%  xent  -  the  horizontal  length  of  the  level  entrance  region. 

%  xtrail  -  the  horizontal  length  of  the  level  exit  region. 

%  Output  parameters 

%  X  -  a  column  vector  containing  a  discrete  amount  of  evenly  spaced  'x* 
%  values  that  were  used  to  calculate  the  profile  and  axle  path. 

%  y  -  a.  column  vector  that  contains  the  'y'  values  for  the  pothole 
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%  profile.  The  bottom  of  the  pothole  is  defined  as  "zero"  y  value. 

%  ya  -  a  column  vector  containing  the  axle  path  'y'  values. 

%  u  -  a  column  vector  containing  the  height  of  the  bottom  of  the 

%  wheel,  'u'  is  the  input  for  the  simulation  model 

%  Other  parameters  assigned  in  the  program: 

%  deltax  -  the  'x'  intraval  spacing. 

%  xent  -  the  horizontal  length  of  the  level  entrance  region. 

%  xtrail  -  the  horizontal  length  of  the  level  exit  region. 

%  rad  -  the  radius  of  the  wheel  that  is  rolling  through  the  pothole. 

%  Varibles  that  are  calculated  based  on  the  above  assigned 
%  parameters: 

% 

%  rhtetal  &  rtheta2  -  the  thetal  &  theta2  in  radians. 

%  xnear  &  xfar  -  the  x  distance  differance  from  the  top  to  the 

%  bottom  for  the  near  &  far  sides. 

%  xcritl  &  xcrit2  -  the  x  value  where  the  wheel  can  no  longer 
%  follow  a  circular  path  as  it  goes  around  the  near  &  far 

%  comers(near  into  hole,  far  out  of  hole). 

%  xcrl  -  the  X  value  where  the  wheel  hits  the  far  side  or  comer 

%  xcr2  -  the  X  value  from  where  the  wheel  hits  the  the  exit  region 

%  from  the  near  side  or  near  comer. 

%  xchl  -  the  X  value  at  the  bottom  comer  on  the  entrance  side. 

%  xch2  -  the  X  value  at  the  bottom  comer  on  the  exit  side. 

%  xch3  -  the  x  value  at  the  top  comer  on  the  exit  side. 

%  xend  -  the  ending  x  value  of  the  pothole  profile. 

%  ycoml  &  ycom2  -  the  y  values  when  the  axle  is  following  a 
%  circular  path  around  the  near  &  far  comers  respectively 

%********assignment  of  pothole  &  wheel  parameters****************** 
deltax=.01 ;  rad=l  7.0; 

%  ♦**♦*♦***  calculating  other  varibles  used  ♦***♦♦*****♦*♦•******♦♦* 
rthetal=thetal  *piyi80;,  rtheta2=theta2*pi/180; 

xnear=depth/tan(rtheta  1 ) ;,  xfar=farht/tan(rtheta2) ; 

xchl=xnear+xent;,  xch2=xchl+xbot; 

xch3=xch2+xfar„  xend=xch3+xtrail; 

xcritl=xent+rad*sin(rthetal);,  xcrit2=xch3-rad*sin(rtheta2); 

xcrit3=xent+rad*sin(rtheta2);,  xcrit4=xch3-rad*sin(rtheta2); 

%  converting/rounding  x-switch  values  to  integer  increments  of 
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%  deltax;  keeps  the  output  x,  y,&ya  vectors  the 

num=xent/deltax;,  nuni=round(num);, 

num  1  =xch  1  /deltax;,num  1  =round(nuni  1 ) ;, 

nuin2:=xch2/deltax;,num2=round(num2);, 

’num3=xch3/deltax;,num3=round(num3);, 

nuin4=xend/deltax;,num4=round(num4);, 

xd  1 =xcrit  1  -xent;,n  1 =xd  1  /deltax;,n  1  =round(  n  1 ) ;, 

nl*deltax; 

xd2=xch3-xcrit2;,n2=xd2/deltax;,n2=round(n2);, 

n2*deltax; 


same  size 
xent=num*deltax; 
xch  1  =num  1  *deltax; 
xch2=num2*deltax; 
xch3=num3*deltax; 
xend=num4*deltax; 
xcritl=:xent  + 

xcrit2=xch3  - 


POTHOLE  PROFILE  *************************** 


%  ******  creating  x  vectors  for  pothole  regions  ***************** 
xl=0:deltax:xent;,  sizl=size(xl); 

x2=:xent+deltax:deltax:xch  1 

x3=xchl  +deltax:deltax:xch2;,  siz3=size(x3); 

x4=xch2+deltax:deltax:xch3;, 
xS=xch3+deltax:deltax:xend;,  sizS=size(xS); 

x=[xl  x2  x3  x4  x5]; 

%  ♦***♦*  calculating  y  values  for  the  pothole  profile  **************** 

yl=depth*ones(l:sizl(2)); 

y2=!  -tan(rthetal)*x2  +tan(rthetal)*xeot  +  depth; 
y3=zeros(l:siz3(2)); 

y4=  tan(rtheta2)*x4  -  tan(rtheta2)*xch2  ; 
y5=  farht*ones(l:siz5(2)); 
y=[yl  y2  y3  y4  y5]; 


^****«***)ti*4i*****  AXLE  PATH  ************************************** 


%******calculating  the  axle  trace  in  the  entrance  region******** 
xal=0:deltax:xent;,sizal=size(xal); 
yal=(depth+rad)*ones(l:sizal(2)); 
if  xcrit2  <=  xcritl, 

%***checking  if  the  axle  motion  is  only  circular***** 
if  xcrit2  >  xent, 

%  finding  the  x  intercept  of  the  near  &  far  side  circular  motion 
for  xint=xcrit2:deltax:xch3, 

yl=depth  +  sqrt((rad/2)  -  ((xint-xent)/2)); 
y2=farht  +  sqrt((rad.''2)  -  ((xint-xch3).''2)); 
if  yl<=y2;  break, end 


no 


end 


for  xint^xent+deltaxrdeltaxixchB, 

yl=depth  +  sqrt((rad/2)  -  ((xint-xent)/2)); 
y2=farht  +  sqrt((rad/2)  -  ((xint-xch3)/2)); 
if  yl<=y2;  break, end 
end 
end 

%  calculating  the  near  and  far  side  y  values  for  the  axle  path 
xa2=xent-*-deltax:deltax:xint; 
ya2=depth  +  sqrt((rad/2)  -((xa2-xent)/2)); 
xa3=xint+de1tax:deltax:xch3; 
ya3=farht  +  sqrt((rad/2)  -((xa3-xch3)/2)); 
yac=[ya2  ya31  ; 

xa4=xch3+deltax:deltax:xend;,  siza4=size(xa4) ; 

ya4=(farht+rad)*ones(  1  :siza4(2)) ; 
ya=[yal  yac  ya4]; 

%checking  &  calculating  axle  path  when  its  not  just  circular 
else 

%****intializes  a  counter  &  the  vector  yac**** 
n=0;,  nnum=(xch3-xent)/deltax;,  yac=2eros(l,nnum); 

%*****checks  &  determination  of  axle  path**** 
for  xasxent+deltax:deltax:xch3, 
if  xa  <=  xcritl, 

ycoml=depth  +  sqrt((rad/2)  -  ((xa-xent)/2)); 

yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2  +rad/cos(rtheta2) ; 
if  ycoml  >  rad  &  ycoml  >  yfar,  yacl=ycoml; 
elseif  yfar  >  rad,  yacl=yfar; 
else  yacl=rad;  end 
elseif  xa  >=  xcrit2, 

ycom2=farht  +  sqrt((Tad.''2)  -  ((xa-xch3)/2)); 
ynear^-tan(rtheta  1 )  *xa+tan(rtheta  1 )  *xent. . . 
+depth+rad/cos(rthetal); 

if  ycom2  >  rad  &  ycom2  >  ynear,  yacl=ycom2; 
elseif  ynear  >  rad,  yacl=ynear, 
else  yacl=rad;  end 
else 

ynear^-tan(rtheta  1  )*xa+tan(rtheta  1  )*xent... 
+depth+rad/cos(rthetal); 

y  far=tan(rtheta2)  *xa-tan(rtheta2)*xch2 +rad/cos(rtheta2) ; 
if  ynear  >  yfar  &  ynear  >  rad,  yacl=ynear; 
elseif  yfar  >  rad,  yacl=yfar; 

1  1  1 


else  yacl=Tad;  end 


end 

n=n+l;,  yac(l,n)=yacl; 
end 

^**m*****m**m**  jyj|g  y  values  in  the  exit  region  **♦*•♦♦♦**♦♦*♦ 
xa4=xch3+deltax:deltax:xend;,  siza4=size(xa4); 
ya4=(farht+rad)  *ones(  1 :  siza4(2)) ; 
ya=[yal  yac  ya4]; 
end 

%**********  PLOT  OF  POTHOLE  PROFILE  &  AXLE  PATH  *♦**♦**♦***♦♦ 

ax=[0,xend,-50,yplot]; 

y=y-depth;  ya=ya-depth;u=ya-17; 

x=x';y=y';ya=ya';u=u'; 

matrix=rx  y  ya  u]; 

save  road2dat.m  matrix 

axis(ax); 

plot(x,y);  hold  on  ;  plot(x,ya)  ,pause,hold  off 
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APPENDIX  E 


NATURAL  FREQUENCIES  AND  MODE  SHAPES 

Table  E-1  below  contains  the  natural  frequencies  and  mode 
shapes  of  the  HUMMV/Howitzer  system.  They  were  calculated  using 
the  procedures  decsribed  in  Chapter  IV.  In  Table  E-1  the  system's 
natural  frequencies  are  listed  in  the  top  row.  The  positional 
coordinate  values  of  the  mode  shape  which  corresponds  to  a 
particular  natural  frequency  are  listed  in  the  column  below  the 
natural  frequency.  The  natural  frequencies  shown  in  Table  E-1  have 
units  of  hertz.  The  positional  coordinates  have  units  of  inches  for  Zi, 
2(2,Z3,  and  Z4,  and  radians  for  the  rotational  cordinates  61,  62,  and  63. 


TABLE  E-1  NATURAL  FREQUENCIES  AND  MODE  SHAPES 

NalJUal _ Frequencies 


Coordi- 

L2.8..6 

1.870 

1.328 

Zi 

-.0001 

.0783 

-.2392 

Z2 

.0130 

.1293 

.1425 

Z3 

.0140 

.2057 

-.0408 

Z4 

-.1064 

.0168 

.0202 

61 

.0002 

.0003 

.0059 

62 

-.0025 

-.0013 

-.0026 

03 

-.0025 

-.0013 

-.0027 

12.26 

I1.M 

18..5J 

23.80 

-.0212 

1.217 

.0130 

.0342 

.0076 

.0364 

-1.222 

-.1221 

.0018 

-.0221 

.0138 

.0277 

-1.385 

-0.221 

-.0097 

.0341 

.0001 

.0011 

.0005 

.0025 

.0001 

-.0003 

-.0001 

-.0036 

.0003 

-.0006 

-.0011 

.0063 
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APPENDIX  F 


A-FRAME  CALIBRATION  DATA 

The  following  tables  contain  the  load  versus  strain  data  that 
was  recorded  during  the  A-Frame  calibration  procedures  that  are 
described  in  Chapter  V.  The  locations  of  the  strain  gages  is  shown  on 
Figures  5-1  and  5-6.  Tables  F-1  and  F-2  contain  the  data  from  the 
vertical  load  tests  and  Tables  F-3  and  F-4  contain  the  data  from  the 
lateral  load  tests.  The  method  used  to  calculate  stress  values  from 
the  strain  gage  readings  is  described  in  Chapter  V. 
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-FRAME  VERTICAL  LOAD  TEST  1 


■FRAME  VERTICAL  LOAD  TEST  2 


TABLE  F-2  VERTICAL  TEST  #2  A-FRAME  CALIBRATION  DATA 


TABLE  F-2  VERTICAL  TEST  #2  A-FRAME  CALIBRATION  DATA  (CONT) 


FRAME  LATERAL  LOAD  TEST  1 


TABLE  F-3  LATERAL  TEST  #1  A-FRAME  CALIBRATION  DATA 


FRAME  LATERAL  LOAD  TEST  1 


-FRAME  LATERAL  LOAD  TEST  2 


APPENDIX  G 


HOWITZER  ORIENTATION  PICTURES  AND  DIAGRAMS 

This  Appendix  contains  pictures  and  diagrams  that  identify  the 
key  components  of  the  Ml  19  Howitzer  that  were  important  for  this 
study. 


Figure  G-1  Towing  in  Firing  Position 


Figure  G-2  Towing  in  Stowed  Position 


Figure  G-3  A-Frame/T-Bar  Support 
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Figure  G-S  Howitzer  Wheels  and  Suspension 
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